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ABSTRACT 


The  process  of  probing  and  measuring  an  HF  channel  for  use  in 
a stored  channel  simulator  will  prove  useful  only  when  the  investi- 
gator has  some  knowledge  of  the  channel  conditions.  Of  particular 
interest  are  the  channel  Doppler  characteristics  and  multipath 
structure.  The  operations  required  to  provide  a measure  of  these 
parameters  include:  (1)  complex  channel  reconstruction;  (2)  pre- 

filtering to  improve  delay  resolution;  (3)  channel  snapshot  genera- 
tion; and  (4)  Doppler  estimation.  A set  of  flexible  Fortran 
programs  which  meet  these  specifications  are  described  in  detail. 
Software  verification  is  achieved  by  means  of  a program  generated 
single  sideband  HF  test  channel.  In  addition,  programming  changes 
to  the  previously  reported  channel  measurement  and  reproduction 
software  are  documented.  These  result  in  significant  decreases  in 
computation  time. 
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SECTION  1 


INTRODUCTION  AND  REPORT  SUMMARY 


The  results  of  probing  and  measuring  an  HF  channel  can  be 
employed,  ostensibly,  to  test  communications  systems  over  actual 
channels  without  resorting  to  full  field  experiment.  If  tests 
are  conducted  with  the  results  of  many  channel  soundings,  a com- 
prehensive picture  of  total  system  performance  may  evolve.  Critical 
to  this  evaluation,  however,  is  an  awareness  of  the  characteristics 
of  the  measured  channel.  Specifically,  channel  definition  by  means 
of  its  multipath  structure  and  Doppler  width  is  important. 

In  Section  3 of  this  report,  we  describe  a software  system 
capable  of  reducing  the  data  on  a channel  probing  tape  to  a mean- 
ingful set  of  parameter  values.  These,  include  delay  power  density, 
Doppler  power  density,  rma  multipath  spread,  rms  Doppler  spread, 
and  mean  Doppler  shift.  In  addition,  modifications  to  the  channel 
measurement  and  reproduction  software  (Part  1 of  this  report),  which 
result  in  less  execution  time,  are  described  in  Section  2. 

Throughout  this  report,  we  frequently  refer  to  a set  of  pro- 
grams which  are  printed  in  Part  1 of  this  final  report,  dated 
15  May  1974.  Each  of  these  programs  can  be  easily  distinguished 
since  they  appear  in  Appendix  P of  that  report  (Figures  P.l  to  P.22). 
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SECTION  2 


CHANNEL  MEASUREMENT  AND  CHANNEL  PLAYBACK  SOFTWARE  MODIFICATIONS 

In  Part  1 of  this  final  report  (covering  the  period  15  October  1973 
to  14  April  1974),  we  discussed  in  some  detail  the  concept  of  a stored 
channel  simulator  capable  of  measuring  and  reproducing  a linear  time- 
variant  HF  channel  over  a finite  bandwidth.  To  demonstrate  the  stored 
channel  technique,  a set  of  Fortran  programs  were  written  which  covered 
four  major  simulator  topics:  1)  probing  signal  generation;  2)  channel 

measurement;  3)  channel  playback;  and  4)  simulator  verification  by 
means  of  a deterministic  time- varying  test  channel. 

Although  no  specific  fidelity  criteria  for  channel  reproduction 
was  adopted,  empirical  evidence  suggested  that  the  deterioration 
resulting  from  all  self-noise  sources  (dc  crosstalk,  harmonic  cross- 
talk, computational  inaccuracies)  was  more  than  70  dB  below  the  desired 
signal.  Unfortunately,  the  programs  took  prohibitively  long  periods 
of  time  to  process  useful  quantities  of  data  (30  seconds  to  1 minute). 

In  the  following  sections  we  will  briefly  describe  the  modifications 
made  to  the  channel  measurement  and  channel  playback  software  to 
achieve  a significant  reduction  in  execution  time. 

2 . 1 Channel  Measurement  Software 

The  program  which  performs  the  channel  measurement  function  is 
entitled  RECDll.  It  is  identical  in  operation  to  the  program  REC2 
as  described  in  the  first  part  of  this  report.  However,  REC2  requires 
approximately  40  hours  of  CPU  time  to  process  one  minute  of  data  — a 
real  time  expansion  factor  of  2400.  In  contrast,  RECDll  will  massage 
the  same  quantity  of  data  in  3.5  hours  - a real  time  expansion  factor 
of  210.  Although  much  of  this  improvement  can  be  attributed  to  NRL' s 
recently  acquired  PDP-11/45  floating  point  processor,  a 50%  time 
savings  has  been  achieved  through  software  changes. 
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A listing  of  RECDll  appears  in  Figure  A.l.  Observe  that,  at 
the  beginning  of  execution,  the  user  is  asked  to  input  the  name  of 
the  output  data  file  he  wishes  to  create  on  magnetic  tape.  Each 
record  of  this  file  will  contain  one  measured  snapshot  (of  length 
255  samples)  of  the  channel  impulse  response.  As  an  example,  suppose 
the  desired  name  is  CHANNL.DAT.  The  user  simply  types  CHANNL  since 
the  program  always  appends  .DAT.  The  only  other  additional  input 
required  by  the  program  is  the  number  of  snapshots  to  be  processed. 

Two  subroutines  are  called  by  RECDll  which  were  not  accessed 
by  REC2.  Subroutine  STASH  (Figure  A. 2)  is  called  with  the  following 
syntax: 


CALL  STASH ( I JKM , ICSDR) 

STASH  simply  places  the  value  of  IJKM  in  the  address  specified  by 
the  variable  ICSDR.  In  this  program,  ICSDR  ■ "177570  which  corres- 
ponds to  the  address  of  the  display  register  on  the  PDP-11/45,  and 
IJKM  contains  the  number  of  snapshots  which  have  been  written  on 
magnetic  tape.  This  gives  the  user  an  indication  of  the  amount  of 
data  processed  without  having  to  resort  to  the  use  of  a device  such 
as  the  keyboard  (whose  device  driver  alone  is  600  words). 

Assembly  language  subroutine  LOOP  (Figure  A. 4)  only  replaces 
seven  lines  of  code  in  REC2.  However,  since  much  of  the  program's 
time  is  spent  executing  those  lines,  their  replacement  by  LOOP  de- 
creased execution  time  by  a factor  of  2.  It  achieves  its  speed  by 
making  very  efficient  use  of  the  floating  point  processor. 

TESTD.DAT  is  the  name  of  the  input  data  file  for  RECDll.  Each 

* 

record  of  the  file  must  contain  510  samples  of  channel  data  (510  sam- 
ples have  a time  duration  equivalent  to  two  times  the  inverse  of  the 
snapshot  rate).  However,  data  derived  from  NRL  channel  probing 
experiments  is  recorded  in  1000-word  12-bit  packed  blocks  on  magnetic 
tape  by  MTIN  (Figure  A. 5).  In  order  to  be  used,  this  data  must  be 
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reformatted  so  that  it  becomes  compatible  with  the  requirements  of 
RECDll.  FORMlA  (Figure  A. 6)  accomplishes  this  objective.  In  addition, 
it  positions  the  data  before  executing  by  skipping  a user  specified 
number  of  files  and  1000-word  records.  Subsequently,  F0RM1A  performs 
the  following  functions: 

(1)  Reads  channel  data  in  1000-word  12-bit  packed  blocks 
from  magnetic  tape  using  call  to  the  NRL  supplied 
input  routine  MTIN.  The  data  is  lef t-justified  in 

a 1000-word  integer  array  IAR. 

(2)  The  data  in  each  wor4  of  IAR  is  shifted  right  by  four, 
bits  for  normalization.  Since  the  data  is  recorded  in 
offset  binary,  "4000  is  subtracted  from  each  word. 

(3)  IAR  is  then  converted  to  floating  point  and  appro- 
priately placed  in  array  X (of  dimension  greater 
than  1000). 

(4)  X is  reformatted  into  510  word  blocks  and  written  on 
magnetic  tape  under  the  file  name  TESTD.DAT. 

FORMlA  replaces  F0RM1  which  appeared  in  Figure  P.5  of  the 
first  report.  As  its  predecessor,  FORMlA  used  modulus  arithmetic 
and  treats  the  data  storage  vector  X as  though  it  were  circular. 

In  this  way,  no  internal  data  moves  need  be  performed  — allowing 
the  programs  to  execute  more  quickly  than  might  normally  be  expected. 
One  additional  feature  has  been  included  which  forces  FORMlA  to  print 
the  magnetic  tape  status  word  on  device  6 should  an  error  occur 
during  a read. 

Function  LSH  (Figure  A. 3)  is  used  by  FORMlA  to  perform  shifting 
operations.  Its  call  must  have  the  syntax  LSH(IW0RD, J) , where  IWORD 
is  an  integer  variable  to  be  shifted  left  by  J bits.  As  an  example, 
for  J“-4,  the  contents  of  IWORD  are  shifted  right  by  four  bits. 
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2 . 2 Channel  Playback  Software 


The  task  of  decreasing  execution  time  for  the  channel  playback 
software  requires  a more  sophisticated  approach  than  that  required 
for  the  channel  measurement  software.  The  approaches  necessarily 
differ  because  the  sources  of  wasted  time  in  the  two  sets  of  software 
differ.  In  the  channel  measurement  case,  certain  Fortran  statements 
were  consuming  huge  amounts  of  CPU  time;  in  the  channel  playback  situ- 
ation, large  amounts  of  time  are  spent  preparing  and  transferring  data 
from  program  storage  to  intermediate  data  storage  on  disk  files.  This 
was  originally  necessitated  by  the  insufficient  core  storage  available 
(20K  words)  on  the  NRL  PDP-11/45. 

The  original  playback  program  PLAY2  is  listed  in  Figure  P.ll  of 
the  first  report.  It  suffers  from  the  extraordinary  complexity 
rendered  by  the  use  of  five  data  files.  This  is  abundantly  obvious 
when  PLAY 2 is  compared  to  a completely  equivalent  program,  PLAY, 
designed  to  run  on  the  PDP-10.  PLAY  is  listed  in  Figure  P.22.  It 
seems  clear  that  the  only  way  to  effectively  speed  program  execution 
is  to  allow  most  intermediate  data  to  remain  in  core.  To  this  end, 
we  have  adopted  an  overlay  program  structure. 

In  order  to  eschew  any  proclivity  toward  prolixity,  we  will 
not  discuss  the  numerous  revisions  necessary  before  working  software 
was  completed.  Merely  observe  that  overlaying  by  itself  is  not 
sufficient.  More  surreptitious  means  for  saving  storage  are  required. 

PLAY2  defines  the  files  to  be  used  by  calls  to  the  Fortran 
subroutine  SETFIL.  SETFIL  is  used  instead  of  ASSIGN  because  random 
access  files  are  needed.  Each  call  to  SETFIL  requires  about  20  words, 
while  SETFIL  itself  is  approximately  250  words  long.  Thus,  just 
setting  the  five  files  in  FLAY2  requires  350  words  of  PDP-11/45  core. 
In  addition,  we  determined  than  only  three  devices  were  required  in 
the  overlayed  program:  1)  PIAYD.DAT  is  the  file  that  contains  data 
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to  be  played  through  the  channel;  2)  RECDT.DAT  contains  the  trans- 
form of  the  recorded  channel  snapshots;  and  3)  PROCD.DAT  contains 
the  processed  data.  In  addition,  the  device  table  on  the  PDP-11/45 
requires  19  words  of  storage  for  each  logical  unit  number  independent 
of  whether  it  is  used.  Since  there  is  nominally  provision  for  eight 
logical  units,  we  could  conceivably  remove  five  logical  unit  numbers 
for  the  overlayed  programs  and  form  what  we  might  term  a minimum 
device  table. 

Taking  this  discussion  to  its  logical  conclusion,  we  can  place 
the  file  names  directly  in  the  minimum  device  table  and  thereby  eli- 
minate the  need  for  program  SETFIL,  five  calls  to  SETFIL,  five  un- 
necessary  logical  unit  numbers,  and  435  words.  A listing  for  the 
device  table  MINDEV,  which  accomplishes  this  objective,  appears  in 
A. 7.  It  assigns  PLAYD.DAT  to  logical  unit  number  1 or  disk  unit  0, 
and  RECDT.DAT  to  logical  unit  number  2 on  magnetic  tape  unit  0.  The 
third  file,  PROCD.DAT,  is  assigned  to  logical  unit  number  3 and  device 
unit  number  1 on  either  disk  or  magnetic  tape.  The  device  chosen 
depends  on  the  output  of  a conditional  assembly.  The  user  required 
alteration  can  be  gleaned  from  the  program  listings. 

Another  source  of  wasted  core  involves  the  Fortran  I/O  package 
and  its  interface  with  the  Fortran  error  handling  routines.  We 
decided  to  put  all  the  device  reads  and  writes  in  the  smallest  over- 
lay segment  and  to  eliminate  any  I/O  from  the  root  segment  and  the 
remaining  overlay  segments.  Unfortunately,  the  Fortran  error  pro- 
cessor requires  a good  deal  of  the  Fortran  I/O  package  and,  of  course, 
it  necessarily  must  appear  in  the  root  segment.  Program  ERRF 
(Figure  A. 8)  allows  us  to  circumvent  this  problem  by  replacing 
the  Fortran  entry  points  ERR,  ERRA,  ERRB,  ERRC,  ERRF  with  a set  of 
entry  points  in  ERRF.  In  case  an  error  occurs,  the  program  puts  a 
message  of  the  form  F030  ABC  XYZ  on  the  keyboard.  In  this  case, 

ABC  ■ error  class,  XYZ  ■ error  number  (both  in  octal).  It  writes  on 
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the  keyboard  without  requiring  a buffer  allocation  by  taking  advantage 
of  a resident  monitor  program.  In  addition,  ERRF  has  another  entry 
point  ERR000  which  can  be  employed  for  additional  error  indications. 
ERFF  saves  over  1000  words  of  storage. 

Figure  2.1  depicts  the  overlay  structure  used  by  the  program. 

The  root  segment  PLAYll  (Figure  A. 9)  has  three  branches.  Branch  B 
contains  subroutine  LAGINT  in  file  LAGl.  LAGINT  appars  in  Figure  P.8 
of  the  first  report.  Branch  C contains  subroutines  SUBl  (Figure  A. 10), 
Z WRITE  (Figure  A. 11),  and  SUB2  (Figure  A. 12).  Branch  D contains  sub- 
routine FFTOVR  (Figure  A. 13)  which,  in  turn,  calls  two  overlay  segments 
Branch  E contains  subroutine  FFT  (Figure  A. 14)  and  Branch  F contains 
subroutine  REALTR  (Figure  A. 15).  A Fortran  equivalent  of  assembly 
language  subroutine  SUB 2 appears  in  Figure  A. 16.  All  of  these  pro- 
grams must  be  linked  together  with  the  aid  of  the  overlay  description 
file  PLAYll. ODL  listed  in  Figure  A. 17. 

As  a result  of  the  programming  changes,  the  channel  playback 
software  now  requires  that  the  data  to  be  played  through  the  channel 
be  formatted  in  255-word  records  instead  of  the  765-word  records 
required  in  the  old  software.  However,  this  requires  only  a change 
in  a parameter  in  program  F0KM2  (Figure  P.12). 
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SECTION  3 


DELAY  SPREAD  AND  DOPPLER  WIDTH  CHARACTERIZATION 
OF  CHANNEL  MEASUREMENT  RECORDINGS 

During  the  course  of  this  program,  two  separate  probings  of  the 
HF  channel  have  resulted  in  several  magnetic  tapes  of  channel  soundings. 
It  is  difficult  to  use  these  tapes  for  other  experimental  work  (such 
as  modem  development)  without  some  knowledge  of  the  character  of  the 
channel  they  represent.  Exclusive  of  additive  noise  interference, 
two  commonly  used  characterization  parameters  are  multipath  spread 
and  Doppler  width.  We  shall  discuss  a set  of  Fortran  programs  which 
are  capable  of  extracting  some  measure  of  these  parameters  from  the 
channel  probing  tapes. 

3.1  System  Background 

Part  1 of  this  report  contains  a discussion  of  the  constraints 
placed  on  the  channel  measurement  process  by  both  the  details  of  the 
NRL  single-sideband  (SSB)  transceiver  and  the  specifics  of  the  probing 
signal  (Sections  2 and  5).  However,  we  can  summarize  their  combined 
effect  on  an  ideally  white  probing  of  the  channel  by  the  spectral 
shape  presented  in  Figure  3.1.  Although  the  long  flat  regions  of 
the  spectrum  allow  channel  characterization  over  a broad  bandwidth, 
the  fast  spectral  rolloff  at  the  band  edges  adds  several  milliseconds 
of  multipath  spread  to  the  multipath  actually  contributed  by  the 
channel . 

For  the  purposes  of  the  channel  measurement  and  playback 
recordings  previously  described,  it  is  important  to  maintain  a 

* 

broadband  "total"  channel  characterization  in  order  that  the  data 

% 

used  during  playback  be  a fair  example  of  the  information  trans- 
mission effects  to  be  encountered.  We  are  interested  in  categorizing 
the  transmission  medium;  therefore,  we  must  reduce  as  much  as  possible 
the  effects  of  the  transceiver  filters  and  the  spectral  weighting  of 
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the  probe  signal.  One  spectral  shape  which  accomplishes  this 
objective  and  which  can  be  applied  to  a baseband  signal,  z(t), 
appears  in  Figure  3.2.  It  is  straightforward  to  show  that  the 
impulse  response  of  this  filter  is 


Xj^t)  - 


sin  ffw^t  cos  ffw^t 
nt  (l  - w^  t^) 


which,  when  sampled  at  a 2w^  rate,  becomes 


(3.1) 


1 ; k-0 
; k-±2 
0 ; elsewhere 


(3.2) 


The  only  additional  processing  required  before  the  channel 
measurement  procedure  can  begin  involves  reconstruction  of  the 
complex  SSB  demodulated  output,  x(t),  where 


z(t ) « Re  x(t)  e 


x! 


J(wo££t+9)' 


(3.3) 


where  0 and  w^^  are,  respectively,  phase  and  frequency  offsets 
between  the  receiver  and  transmitter  filters.  Assuming  that  w^^ 
is  small  compared  to  any  dead  zone  around  dc  in  Figures  3.1  and  3.2, 
then  we  can  say  that 


J(w  ff+8) 

z(t)  + j z(t)  • x(t)  e 


(3.4) 


where  z(t)  is  the  Hilbert  transform  of  z(t). 

Of  course,  we  could  have  formed  the  complex  signal  z(t)+j  z(t) 
before  weighting  by  the  SSB  equivalent  of  the  spectrum  in  Figure  3.2. 
However,  software  considerations  were  preeminent  in  the  rejection 
of  this  procedure. 
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3.2  Software 


A block  diagram  of  the  required  software  is  given  in  Figure  3.3. 
Observe  that  we  have  decreased  the  order  of  the  Interpolator  used 
during  channel  measurement  from  a fourth-order  Lagrange  to  a second- 
order  Lagrange.  This  was  necessitated  by  the  interplay  between  the 
higher  core  requirements  for  complex  processing  and  the  finite 
amount  of  core  available.  However,  since  the  requirements  for 
channel  characterization  are  less  stringent  than  those  for  channel 
reproduction,  the  interpolator  order  reduction  is  not  significant. 

For  example,  in  Part  1 of  this  report  we  noted  that  0.1  dB  Doppler 
attenuation  occurs  at  +5.65  Hz  for  a fourth-order  interpolator  and 
at  +2.4  Hz  for  a second-order  interpolator.  At  5.65  Hz,  a second- 
order  interpolator  renders  only  0.6-dB  attenuation  — an  acceptable 
quantity  for  channel  characterization. 

The  Hilbert  transformer  depicted  in  Figure  3.3  can  be  imple- 
mented in  a great  many  ways.  We  chose  to  use  a finite  impulse 
response  (FIR)  linear  phase  digital  filter.  Toward  this  end,  we 
employed  a program  written  by  McClellan,  Parks,  and  Rabiner  for 
optimum  FIR  filters  [3.1].  The  design  was  restricted  by  three 
criteria:  1)  that  the  frequency  response  of  the  filter  appear 

white  with  respect  to  the  transceiver  filter-probing  signal  spec- 
tral combination  of  Figure  3.1;  2)  that  passband  ripple  magnitude 
be  less  than  1%;  and  3)  that  the  order  of  the  filter  (the  duration 
of  its  impulse  response)  be  odd.  (Since  for  N odd  every  other 
sample  is  zero,  the  last  requirement  is  important  to  minimize 
computation.)  A FIR  filter  of  length  43  satisfies  these  criteria. 

The  Hilbert  transformer  is,  of  course,  odd  symmetric  and  unre- 
alizable. As  a result  of  its  finite  Impulse  response,  the  filter 
can  be  made  realizable  by  inserting  a delay  of  21  samples  in  the 
untransformed  path.  This  is  Indicated  in  Figure  3.3.  The  odd 
symmetry  is  computationally  important  since  it  allows  the 


programmer  to  replace  two  multiplies  by  a multiply  and  a subtract. 
Thus,  the  43  sample  impulse  response  requires  only  11  multiplications 
and  11  subtractions. 

Figure  3.4  is  a block  diagram  of  a technique  for  measuring  the 
spectral  characteristics  of  a Hilbert  transformer.  The  amplitude 
response  of  the  filter  on  linear  and  log  scales  appears  in  Figures  3.5 
and  3.6,  respectively.  Observe  the  equiripple  response  of  this 
optimum  filter  in  both  the  passband  and  stopband. 

We  chose  to  combine  (1)  prefiltering  to  improve  delay  resolution, 
(2)  complex  channel  reproduction,  and  (3)  the  channel  measurement 
procedure  into  one  program  named  RECSCF  (A.  18).  The  word  RECSCF  is 
a slightly  misnamed  acronym  for  RECord  Scattering  Function.  (A 
scattering  function  is  not  actually  calculated  — only  an  estimate 
of  its  transform  along  the  Doppler  dimension.) 

Program  RECSCF  requires  four  files:  1)  a random  access  file 

for  temporary  data  storage  STORl.DAT;  2)  a file  LAGR.DAT  containing 
second-order  Lagrange  Interpolator  coefficients  [created  by  SETLAG 
(A. 19)];  3)  a file  TESTD.DAT  formatted  in  records  of  510  words  con- 
taining the  data  to  undergo  the  channel  measurement  process;  and 
4)  a magnetic  tape  output  file  whose  name  is  inputted  by  the  user. 

In  addition,  the  user  must  specify  the  number  of  snapshots  to  be 
processed  by  the  RECSCF.  If  the  number  of  snapshots  specified  is  • 
greater  than  the  length  of  the  input  data  file,  the  last  record 
read  is  processed,  a normal  termination  ensues,  and  the  user  is 
informed  that  an  end  of  file  has  been  encountered.  The  program 
performs  one  subtle  maneuver  to  save  execution  time  by  subtracting 
"200  from  the  first  word  of  a floating  point  number  instead  of 
dividing  by  2. 

When  the  program  terminates,  a summary  of  the  total  power  at 
each  tap  is  written  on  the  line  printer.  This  information  is  im- 
portant to  the  user  who  must  subsequently  determine  which  taps  are 
to  be  examined  for  Doppler  spread  characteristics. 
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SAMPLE 


HILBERT 

TRANSFORM 


FFT 


Figure  3.4  Hilbert  Transformer  Test 


Define  g(kT,nA)  to  be  one  complex  snapshot  of  the  channel 
impulse  response,  where  ;jr  is  the  snapshot  rate,  and  ~ is  the 
sampling  rate  (2w^  in  Figure  3.2).  As  discussed  in  Part  1 of  this 
report,  each  tap  in  a tap  delay  line  channel  model  is  separated  by 
A seconds.  During  playback,  each  tap  coefficient  is  updated  every 
T seconds.  As  an  indication  of  delay  spread,  we  are  interested  in 
the  power  in  each  tap  averaged  over  the  number  of  snapshots  processed. 
Normally,  only  a few  taps  will  contain  significant  energy.  For 
example,  a 1-ms  multipath  spread  at  a 9-kHz  sampling  rate  would 
involve  only  9 taps.  However,  coupled  with  a bandwidth  reduction 
filter  length  of  five  samples  (3.2),  this  yields  13  taps  which  con- 
tain some  meaningful  information.  The  starting  and  end  indices  of 
the  taps  with  significant  energy  can  be  ascertained  from  the  RECSCF 
line  printer  output.  It  is  on  this  set  of  taps  that  we  would  like 
to  perform  spectral  estimation  in  the  Doppler  dimension. 

Techniques  for  calculating  power  spectra  are  straightforward 
and  well  cataloged.  The  one  which  we  selected  for  this  software 
system  falls  under  the  general  topic  of  linearly  modified  spectral 
estimates  [3.3].  This  method  is  distinguished  by  the  subdivision 
of  data  into  windowed  segments  prior  to  transforming.  A succinct 
description  follows  of  the  steps  involved  in  producing  that  which 
might  be  termed  in  the  ideal  case  "consummate  spectra". 

(1)  Choose  the  length  of  the  FFT  to  be  used  according  to 
the  Doppler  resolution  required.  For  example,  fre- 
quency resolution  of  an  FFT  “ “ Hz,  where  N is  the 
number  of  samples  in  the  FFT  and  T is  the  distance 

between  time  samples  of  the  FFT  data  base.  In  this 
% 1 

case,  ~ (35.294  Hz)  corresponds  to  the  inverse  of 
the  snapshot  rate.  For  N-64,  frequency  resolution 
is  approximately  0.551  Hz,  which  for  many  applications 
is  adequate.  Of  course,  we  could  increase  the  length 
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of  the  FFT  to  attain  better  resolution  if  we  are 
willing  to  increase  the  variance  of  the  spectral 
estimate  [3.4]. 

(2)  Choose  a time  window  for  weighting  the  data.  There 
are  obviously  many  from  which  to  choose.  In  addition, 
once  the  weights  have  been  established,  all  are 
equally  easy  to  apply.  The  one  we  have  chosen  to 

use  is  termed  an  optimal  window  by  A.  Eberhard  [3.2]. 
This  window  maximizes  the  ratio  of  the  energy  in  the 
mainlobe  of  the  window  spectrum  to  the  total  energy 
in  the  spectrum. 

(3)  Apply  overlapping  windows  to  the  data  base,  as  indi- 
cated in  Figure  3.7.  This  method  minimizes  the  loss 
in  potential  spectral  information  that  nonequal 
weighting  of  the  data  base  renders.  Ideally,  the 
sum  of  the  values  of  all  windows  in  the  data  base 

at  each  element  in  the  data  base  should  be  unity  [3.5]. 

(4)  Weight  each  overlapping  segment  in  the  data  base  with 
the  window  coefficients,  take  a FFT,  find  the  magni- 
tude square,  and  average  over  all  segments.  For  the 
parameter  values  listed  above,  there  are  approxi- 
mately 45  periodograms  averaged  for  each  tap,  given 
one  minute  of  input  data. 


A program  for  computing  Doppler  spectra  at  each  tap  position, 
DDSPEC,  is  listed  in  Figure  A. 20.  The  user  must  input  the  name  of 
the  magnetic  tape  data  file  which  contains  the  recorded  snapshots 
of  the  channel  impulse  response.  He  must  also  enter  the  indices 
of  the  left  and  right  tap  boundaries  surrounding  the  taps  to  be 
processed.  For  example,  if  the  user  wants  taps  245-255  and  1-6 
processed  (17  taps),  he  merely  enters  245,6,.  The  maximum  circular 
separation  of  the  tap  boundaries  entered  must  be  less  than  or  equal 
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to  84.  In  the  example  just  given,  the  circular  tap  separation  is 
16.  The  program  creates  a random  access  output  data  file  DELDOP.DAT 
which  contains  the  final  spectral  estimates  at  each  tap  position 
specified.  A summary  of  this  file  is  written  on  the  line  printer 
just  prior  to  program  termination. 

Another  program,  DDREAD  (Figure  A. 21),  massages  DELDOP.DAT  and 
prints  its  contents  in  a format  which  is  conducive  to  viuual  extrac- 
tion of  gross  physical  details  of  the  delay  Doppler  surface.  Up  to 
14  taps  can  be  presented  in  this  way  at  any  one  time.  An  example  of 
this  format  will  be  given  in  the  next  section. 

For  the  user  who  would  rather  see  a statistical  summary  of  the 
information  collected,  program  DDSTAT  (Figure  A. 22)  is  available. 

It  reports  the  following  information  via  the  line  printer: 

1)  marginal  delay  power  distribution  for  the  taps  examined  and 
the  time  in  seconds  by  which  each  tap  is  separated;  2)  marginal 
Doppler  power  distribution  over  a specified  range  of  frequencies; 

3)  mean  delay  with  respect  to  the  first  tap  examined;  4)  rms  delay; 

5)  mean  Doppler;  6)  rms  Doppler  spread;  7)  the  taps  at  which  delay 
peaks  occur  and  time  differences  In  seconds  between  peaks;  and 
quencies  in  Hz  at  which  Doppler  peaks  occur;  9)  Doppler  statistics 
at  each  tap  delay;  and  10)  graphs  of  marginal  delay  power  and  marginal 
Doppler  power.  An  example  Is  given  in  Section  3.4. 

3.3  Software  Verification  (Basic  Programs) 

For  the  purposes  of  software  verification,  it  is  necessary  to 
create  a complex  SSB  HF  channel  with  typical  multipath  and  Doppler 
characteristics.  The  approach  we  shall  take  is  summarized  in 
Figure  3.8.  This  system  is  designed  to  emulate  probing  with  a PN 
sequence  — the  result  of  which  is  a set  of  frequency  samples  of  the 
HF  channel  where  the  intersample  space  in  Hz  is  equal  to  the  Inverse 
of  the  length  in  seconds  on  one  period  of  the  probing  signal.  (More 
details  of  the  probing  process  are  presented  in  Part  1,  Section  3, 
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of  this  report.)  An  equivalent  software  approach  involves  weighting 
the  line  spectrum  of  a periodically  repeated  PN  sequence  by  a simu- 
lated complex  HF  channel. 

Ideally,  the  response  of  one  path  of  the  HF  channel  to  a 
sinusoidal  input  of  the  form 

x(t)  ■ <s/2P  sin(wt  + 0)  (3.5) 


is 


y(t)  " JT?  r cos(wt  + g +<p) 


(3.6) 


where  r is  the  time-varying  Rayleigh-distributed  amplitude  of  the 
path  and  <p  is  the  phase  resulting  from  Doppler, time  delay,  etc.  In 
general,  a complex  noise  process  should  be  added  to  y ( t ) ; however, 
for  the  purposes  of  test  channel  generation,  we  shall  only  consider 
the  noiseless  case, 

Similarly,  the  received  signal  over  a multipath  channel  at  any 
tone  is 


y(t) 


r^  cos  (wt  + 0 +<Pj) 


(3.7) 


where  the  r^  are  independent  Rayleigh-distributed  variables  and  N 
is  the  number  of  paths.  For  normalization  purposes,  we  will  con- 
strain  the  sum  over  the  power  in  the  j path  (Pj)  to  be  unity: 


(3.8) 
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For  an  array  of  M tones 


y(t) 


cos^t  +ei  +0iJ) 


(3.9) 


This  model  assumes  that  is  constant  over  all  frequencies  of 
interest  while  <p^j  is  frequency-dependent.  Both  r^  and  the  Doppler- 
dependent  portion  of  can  be  obtained  from  a set  of  filtered 
complex  Gaussian  variables.  In  the  program  to  be  described,  fil- 
tering is  accomplished  with  a two-pole  filter  whose  passband  shape 
approximates  typical  Doppler  rolloff  characteristics. 

A program,  HFSIMS,  which  creates  a SSB  HF  channel  by  frequency 
weighting  a set  of  tones  according  to  the  desired  multipath  and 
Doppler  characteristics  is  listed  in  A. 23.  The  output  magnetic  tape 
data  file,  HFSIMS.DAT,  will  contain  transforms  of  the  baseband  snap- 
shots. Inputs  to  HFSIMS  include:  1)  the  number  of  snapshots  to  be 

written  on  magnetic  tape;  2)  the  number  of  frequency  samples  in  the 
sideband  of  interest  (by  default,  the  program  chooses  the  upper  side- 
band); 3)  the  frequency  separation  between  tones;  4)  the  snapshot 
rate  (because  of  the  details  of  the  probing  process,  items  3 and  4 
must  be  equal);  5)  the  number  of  paths  in  the  channel  model  (non- 
impulsive  path  delay  densities  are  not  computationally  feasible 
since  any  simulation  using  that  model  would  require  huge  amounts  of 
CPU  time);  6)  the  Doppler  bandwidth  of  each  path  (Hz);  7)  the  power 
of  each  path  in  dB  (for  normalized  results,  the  total  power  in  the 
paths  should  be  unity);  8)  the  Doppler  offset  of  each  path  (Hz); 
and  9),  the  relative  delay  (seconds)  of  each  path  (the  first  path 
should  have  a relative  delay  of  0). 
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Armed  with  synthetic  snapshots  of  an  HF  channel,  we  are  ready 
to  simulate  probing  the  channel  with  a periodically  repeated  PN 
sequence.  This  situation  is  identical  to  one  encountered  in  Part  1 
of  this  report  in  which  it  was  necessary  to  play  data  through  re- 
corded snapshots  of  the  channel.  Program  PLAYll  (A. 9)  accomplished 
this  objective.  It  seems  reasonable  then  to  adopt  PLAYll  for 
"playing  back"  the  PN  probing  signal  through  the  synthetic  snapshots. 

One  simplification  is  in  order.  PLAYll  is  capable  of  handling 
large  bandwidth  Doppler  spreads  or  large  Doppler  shifts  without  im- 
parting significant  passband  deterioration  (Section  2 and  3 of 
Part  1).  This  is  a result  of  using  a fourth-order  Lagrange  inter- 
polator for  zeroing  Doppler  images. at  multiples  of  the  snapshot 
rate.  Since  we  are  only  testing  the  software,  it  is  not  necessary 
for  these  initial  experiments  to  use  relatively  large  Dopplers. 

As  a result,  it  will  prove  computationally  expedient  to  use  only  a 
second-order  interpolator  while  limiting  maximum  Doppler  shifts  to 
around  3 Hz.  PLYSCF  incorporates  all  these  features  into  one 
program. 

PLYSCF  assigns  four  different  data  files  during  the  course  of 
its  execution.  File  LAGR.DAT  is  created  by  program  SETLAG  (Figure 
A. 19).  File  PNSEQT.DAT  contains  a 768  point  real  transform  of  one 
period  of  a PN  sequence  of  length  255.  Program  SETPNS  (Figure  A. 25) 
creates  this  file  while  employing  the  use  of  an  FFT  program  speci- 
fically modified  to  perform  a 384  point  complex  transform  [FFT  384 
(Figure  A. 26)].  File  HFSIMS.DAT  contains  the  transformed  channel 
snapshots  generated  by  HFSIMS . The  last  file,  TESTD.DAT  contains  the 
result  of  playing  a PN  sequence  through  the  synthetic  channel. 

Should  an  end  of  file  occur  in  HFSIMS.DAT,  the  program  will  termi- 
nate normally  after  printing  the  record  number  being  processed 
when  the  EOF  occurred.  PLYSCF,  like  most  of  the  main  programs 
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previously  described,  puts  a numerical  indication  of  the  record 
being  processed  on  the  display  register. 

An  overlay  program,  PLYOVR  (Figure  A. 27),  is  functionally  equi- 
valent to  PLYSCF . It  has  the  capability,  however,  of  being  able  to 
incorporate  any  order  interpolator  by  changing  the  value  of  the 
parameter  IQ  to  reflect  the  order  desired.  Some  obvious  changes 
must  also  be  made  in  array  dimensioning.  The  user  can  also  save 
an  additional  400  words  by  incorporating  a variation  of  the  minimum 
device  table  program,  MINDEV.  The  overlay  segments,  SUBl,  FFTOVR, 
and  SUB2  are  listed  in  Figures  A. 28,  A. 29,  and  A. 30,  respectively. 
PLYOVR. ODL,  which  contains  the  necessary  overlay  descriptor  language, 
appears  in  Figure  A. 31. 

PLYOVR  uses  both  synchronous  and  asynchronous  manual  load  opera- 
tions, as  specified  by  the  PDP-11/45  linker  manual.  (It  should  be 
observed  that  the  bit  used  as  a synchronous-asynchronous  switch  in 
CALL  LOAD  is  reversed  in  all  the  documentation.)  Manual  load  saves 
both  time  and  memory  when  compared  to  the  performance  of  the  PDP-11/45 
autoload  feature.  PLYOVR,  using  a second-order  interpolator,  executes 
quite  rapidly i Approximately  two  seconds  of  computation  per  snapshot 
are  required  — corresponding  to  a real-time  expansion  factor  of  70. 

3.4  Test  Examples 

To  establish  program  credibility,  we  have  devised  test  examples 
which  will  exercise  many  of  the  software  features.  A summary  of  the 
programs  required  in  the  order  of  their  execution  appears  in 
Table  3-1.  As  previously  discussed,  SETPNS  and  SETLAG  create 
data  files  necessary  for  the  synthetic  test  channel.  The  program 
which  determines  the  characteristics  of  that  channel  is  HFSIMS. 
Following  the  stated  guidelines,  we  chose  to  keep  the  maximum  ex- 
cursion of  significant  Doppler  energy  to  3 Hz.  Also,  to  insure 
that  Doppler  estimation  would  be  performed  over  many  Doppler  cycles 
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| TABLE  3-1 

SUMMARY  OF  PROGRAMS  REQUIRED  FOR  TEST  EXAMPLE 

\ 


Program 

Input  Files 

Output  Files 

Comments 

SETPNS 

PNSERT.DAT 
(Random  Access) 
(DK) 

Put  transform  of  3 repe- 
titions of  PN  sequence 
in  file  PNSEQT.DAT. 

SETLAG 

m m 

LAGR.DAT 

(DK) 

LAGD.DAT 

(DK) 

Lagrange  interpolation 
coefficients 
LAGR.DAT  - order  2 
LAGD.DAT  - order  4 ■ 

HFSIMS 

am  m 

HFSIMS.DAT 

(DK) 

Put  transform  of  channel 
snapshots  in  HFSIMS.DAT. 

PLYSCF 

LAGR.DAT 

(DK) 

PNSEQT.DAT 
(Random  Access) 
(DK) 

HFSIMS . DAT 
(MT) 

TESTD.DAT 

(DK) 

i 

Create  test  channel 
(TESTD.DAT)  by  playing 
PN  sequence  through 
channel  snapshots 
created  by  HF  simulator 

RECSCF 

LAGR.DAT 

(DK) 

TESTD.DAT 

(DK) 

STORl . DAT 
(DK) 

NAME . DAT 
(MT) 

TEST.DAT  contains  base- 
band signal.  Program 
prefilters,  creates  com- 
plex channel,  correlates, 
and  interpolates  with 
Lagrange  interpolator  of 
order  2.  Records  of 
NAME. DAT- 25 5. 

CHFFT 

NAME . DAT 
(DK) 

FFT255.DAT 

(DK) 

Read  in  record  M of 
NAME . DAT  FFT  and  put  in 
file  FFT255.DAT. 

FFT255.DAT 

(DK) 

-- 

Plots . 

DDSPEC 

NAME . DAT 
(MT) 

DELDOP . DAT 
(DK) 

Compute  delay  Doppler 
spectrum. 

even  at  the  lowest  Doppler  resolution  Increment  (0.551  Hz),  we 
chose  to  perform  the  first  test  over  40  seconds  of  data.  -As  a 
first-order  estimate,  this  implies  that  spectral  averaging  will 
include  at  least  22  independent  samples.  Therefore,  we  can  expect 
the  standard  deviation  of  any  statistical  results  to  be  21.37.  of 
the  actual  value. 

User  input  to  HFSIMS  included:  1)  CHANNL  as  the  name  of  the 

output  data  file  (.DAT  assumed);  2)  1410  (40  second)  snapshots  to 
be  processed;  3)  128  frequency  samples  in  the  upper  sideband 
(including  one  sample  at  the  lower  band-upper  band  transition 
point);  4)  a 35.294-Hz  snapshot  rate;  and  5)  a separation  between 
frequency  samples  of  35.294  Hz.  Data  specified  for  a two-path 
test  is  presented  in  Table  3-2. 

After  executing  HFSIMS,  FLYSCF,  and  RECSCF,  the  reconstructed 
synthetic  channel  snapshots  are  written  in  file  CHANNL.DAT.  As  an 
intermediate  check  on  software  execution,  programs  CHFFT  (Figure  A. 32) 
and  P1XXTCH  (Figure  A. 33)  plot  the  magnitude  in  dB  of  an  FFT  of  any 
desired  record  in  CHANNL.DAT.  NRL's  Versatec  matrix  plotter  soft- 

4.1. 

ware  is  required  for  PLOTCH.  A plot  of  the  14*"  such  record  appears 
in  Figure  3.9.  Although  the  exact  features  of  the  two-path  model 
are  not  easily  discerned  because  they  have  been  all  but  obliterated 
by  the  effects  of  Doppler,  some  important  observations  can  be  made. 

On  the  frequency  normalized  to  sampling  rate  scale,  the  upper  side- 
band falls  in  the  region  0.0  to  0.5.  Notice  the  small  amount  of 
energy  in  the  lower  sideband.  Both  sidebands  were  prefiltered  in 
RECSCF  with  a raised  cosine  weighting.  The  nulls  of  this  weighting 
at  normalized  frequencies  0.  and  .5  are  obvious.  In  addition,  the 
raised  cosine  weighting  is  apparent  in  the  mainlobe  of  the  upper 
sideband. 


TABLE  3-2 


ir 


L 


TWO- PATH  TEST  DATA  SUMMARY 


•i 

j 

$ 

A 


Parameter 

1 

Path  1 

1 

Path  2 

Doppler  Standard 
Deviation  (Uz) 

0.5 

1.5 

(RMS  Doppler  Spread) 

(1.0) 

(3.0) 

Doppler  Shift  (Hz) 

1.5 

-1.5 

Path  Delay  (sec) 

0 

0.0005 

9* 


In  order  to  more  graphically  illustrate  the  influence  of  a 
two-path  model  on  spectral  particulars,  we  performed  a separate 
test  in  which  neither  path  was  subject  to  any  Doppler  effects.  The 
results  of  plotting  the  spectrum  of  the  14  record  of  CHANNL.DAT 
appears  in  Figure  3.10.  The  characteristics  to  which  we  referred  in 
the  previous  paragraph  are  vividly  portrayed. 

DDSPEC  completes  the  basic  software  processing.  In  addition 
to  the  name  of  the  input  data  file,  it  requires  the  indices  of  the 
left  and  right  tap  boundaries  of  the  taps  to  be  examined.  This  can 
be  obtained  from  RECSCF’s  printout  of  power  at  each  tap  location 
summed  over  all  snapshots.  The  data  for  this  test  is  reproduced  in 
Figure  3.11.  Observe  that  nearly  all  delay  power  falls  between 
taps  20  and  33.  After  running  DDSPEC,  DDREAD  should  be  executed 
(if  necessary,  several  times  until  all  the  taps  processed  are  ex- 
amined). Because  of  space,  its  output  is  given  in  two  parts. 

Figure  3.12(a)  contains  Dopplers  from  -17.647  Hz  to  -0.551  Hz; 

Figure  3.12(b)  contains  Dopplers  from  0 Hz  to  17.096  Hz.  The 
peaks  of  the  delay  Doppler  lobes  for  each  path  are  underlined. 
Observe  that  they  are  separated  by  4-1/2  taps  (0.0005  second). 

The  user  can  get  a feeling  for  spread  in  the  Doppler  dimension  by 
comparing  the  ratio  of  Doppler  power  densities  at  two  similar  points 
for  each  path.  In  this  manner,  it  is  easy  to  see  that  the  second 
path  has  a much  broader  Doppler  bandwidth  than  the  first  path. 
However,  the  peak  of  the  second  path  is  smaller  than  that  of  the 
first  because  the  integrated  power  density  for  the  two  paths  was 
user-specified  to  be  equal. 

A statistical  summary  of  the  data  in  Figure  3.12  is  given  by 
the  output  of  DDSTAT  in  Figure  3.13.  Observe  that  the  claims  for 
delay  peak  and  Doppler  peak  locations  are  proportional  to  some 
multiple  of  an  index  of  delay  or  Doppler  increments,  For  example, 
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FREQUENCY  NORMALIZED  TO  SAMPLING  RATE 


t~h 

Figure  3.10  FFT  Magnitude  of  14  Snapshot  of  RECSCF  Output 
(Two-Path  Model  without  Doppler) 
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Figure  3.11  Tap  Power  Summed  Over  1410  Snapshots • for  First  Test  Example 
(RECSCF  Output) 
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DOPPLER  STATISTICS  FOR  THE  14  TAPS  EXAMINED: 
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Figure  3.13  (Continued) 


a Doppler  peak  at  -1.654  Hz  corresponds  to  three  times  the  Doppler 
step  size  of  0.551  Hz,  while  the  actual  Doppler  peak  is  at  -1.5  Hz. 
More  accurate  peak  location  by  numerical  interpolation  can  be  in- 
cluded if  the  need  arises  for  more  than. just  gross  channel  parameters. 

Of  particular  interest  in  Figure  3.13  are  the  Doppler  statistics 
for  the  14  taps  examined.  According  to  Table  3-2,  the  channel  path 
at  0 relative  delay  should  have  an  rras  Doppler  spread  of  1 Hz,  while 
the  path  at  0.5  msec  relative  delay  should  have  an  rms  Doppler 
spread  of  3.0  Hz.  Noting  from  the  printout  that  delay  peak.s  occur 
at  taps  5 and  9,  we  can  easily  verify  that  the  measured  rms  Doppler 
spreads  for  the  two-paths  are  1.28  Hz  and  2.93  Hz,  respectively. 

These  numbers  seem  reasonable  in  view  of  the  previously  discussed 
21%  standard  deviation  in  measurement. 

» 

As  a visual  aid  to  the  user,  DDSTAT  plots  marginal  delay  power 
and  marginal  Doppler  power.  These  appear  in  Figures  3.14  and  3.15, 
respectively.  Observe  that  the  0.5-msec  path  separation  is  easily 
resolved  by  the  programs.  The  3-Hz  Doppler  shift  between  paths  is 
also  easily  resolved.  In  fact,  it  seems  reasonable  that,  for  small 
Doppler  spreads,  less  than  1-Hz  resolution  should  be  possible  for 
Doppler  peaks  at  the  same  delay. 

In  an  effort  to  more  clearly  demonstrate  the  resolution  capabili- 
ties of  the  software  processor,  we  have  conducted  a second  test  with 
each  channel  path  less  distinct  from  the  others  in  delay  and  Doppler 
than  the  channel  paths  in  the  first  test.  A summary  of  the  delay 
and  Doppler  characteristics  for  this  three-path  test  case  is  pre- 
sented in  Table  3-3.  The,  relative  path  delays  specified  are  chosen 
to  place  peaks  in  the  measured  delay  power  spectrum  at  integral 
multiplies  of  the  intertap  spacing  in  seconds.  This  allows  best- 
case  estimation  of  software  delay  resolution  capabilities.  Observe 
that  the  intertap  spacings  correspond  to  2 and  3 taps,  respectively  — 
in  distinction  to  the  first  test  in  which  the  path  delay  was  equiva- 
lent to  4-1/2  taps. 
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TOP  NUMBER 


Figure  3.14  Marginal  Delay  Power,  in  dB  for  First  Test  Example 
(Each  tap  separated  by  0.111  msec) 
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Figure  3.15  Marginal  Doppler  Power  for  First  Test  Example 
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[ TABLE  3-3 

| THREE -PATH  TEST  DATA  SUMMARY 

; I 

: i ■ : 


Parameter 

Path  1 

Path  2 

Path  3 

Doppler  Standard 
Deviation  (Hz) 

■a 

mm 

0.25 

(RMS  Doppler  Spread) 

(0.5) 

Power  (dB) 

-4.771 

-4.771 

-4.771 

Doppler  Shift  (Hz) 

-1.102 

-1.653 

Path  Delay  (sec) 

' 0 

0.000222 

0.000555 

Doppler  shift  magnitudes  are  constrained  to  be  (1)  less  than 
3 Hz  in  order  to  minimize  deterioration  resulting  from  interpola- 
tion  filters  in  RECSCF  and  PLYSCF,  and  (2)  greater  than  1 Hz  in 
order  that  spectral  estimates  be  averaged  over  as  large  a number 
of  independent  samples  as  possible.  Since  the  test  has  been  con- 
ducted over  60  seconds  of  data,  even  at  0.55  Hz  we  can  expect  a 
standard  deviation  in  measurement  of  177*.  In  addition,  the  mean 
of  each  Doppler  lobe  is  placed  at  an  integer  multiple  of  resolution 
cell  increments  (0.551  Hz).  The  rms  Doppler  spread  of  each  path  is 
less  than  one  resolution  cell. 

As  noted  previously,  RECSCF  (in  addition  to  its  function  of 
measuring  the  complex  channel  snapshots)  also  outputs  the  average 
tap  power  for  the  entire  processing  time.  This  information  is  re- 
produced in  Figure  3.16.  The  two  measured  peaks  are  underlined. 
Without  resolution  limitations,  the  three-path  peaks  would  have 
appeared  at  the  24th,  26th,  and  29th  tap.  It  is  clear  that  the 
peak  which  occurs  at  the  25  tap  is  a result  of  the  first  and 
second  paths  combining. 

The  results  of  running  DDSPEC  and  DDREAD  appear  in  Figures 
3.17(a)  and  3.17(b).  As  before,  delay  Doppler  peaks  are  underlined. 
It  is  relatively  easy  to  determine  peaks  in  the  two-dimensional 
delay  Doppler  surface  before  marginal  distributions  are  derived. 
However,  the  user  must  resort  to  first  derivative  estimation  in 
order  to  distinguish  the  first  two  paths. 

The  statistical  information  obtained  with  the  aid  of  DDSTAT 
is  presented  in  Figure  3.18.  As  expected,  the  program  finds  only 
two  peaks  in  delay  and  Doppler.  It  should  be  remembered,  however, 
that  the  number  of  peaks  in  the  marginal  distributions  is  not 
necessarily  the  same  as  the  number  of  distinguishable  peaks  in 
the  delay  Doppler  surface. 
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Figure  3.17(a)  Doppler  Power  vs.  Delay  Power,  -17.65  to  -0.551  Hz 

(Second  Test  Example) 
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(Second  Test  Example) 
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Figure  3.18  Statistical  Summary  of  Second  Test  Example  (DDSTAT  Output) 
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Figure  3..  18  (Continued) 


The  Doppler  statistics  for  the  14  taps  examined  merit  some 
comment.  The  measured  Doppler  shifts  at  taps  4,  6,  and  10  of 
-1.13  Hz,  -1.4  Hz,  and  -2.71  Hz,  respectively,  agrees  well  with 
the  specified  values  of  -1.1  Hz,  -1.65  Hz,  and  -2.76  Hz,  especi- 
ally considering  the  measurement  corruption  resulting  from  path 
interference.  However,  the  measured  Doppler  spread  at  each  of 
these  taps  is  approximately  twice  the  value  which  had  been  expected. 
This  is  a result  of  the  mainlobe  width  of  the  window  function  used 
for  spectral  estimation.  This  can  be  reduced  by  increasing  the 
length  of  the  window  (and  the  associated  FFT) ; however,  for  a fixed 
duration  data  sample,  this  would  reduce  the  number  of  periodograms 
over  which  averaging  is  performed  and  increase  the  variance  of  the 
estimate. 

Graphical  presentation  of  marginal  delay  power  and  marginal 
Doppler  power  appears  in  Figures  3.19  and  3.20,  respectively. 
Coalescence  of  the  first  two  paths  in  both  delay  and  Doppler  is 
obvious.  However,  the  position  of  the  third  path,  which  is  sepa- 
rated from  the  others  by  a minimum  of  three  taps  and  1.1  Hz,  is 
easily  extracted.  We  might  estimate  the  lower  resolution  of  the 
present  software  system  to  be  ■j  msec  in  delay  and  1 Hz  in  Doppler. 
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TOP  NUMBER 


Figure  3.19  Marginal  Delay  Power  (Second  Test  Example) 


Figure  3.20  Marginal  Doppler  Power  (Second  Test  Example) 
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APPENDIX  A 


PROGRAM  LISTINGS 


DIMENSION  X(255) ,Y<509) ,Z< 1020) ,H< 1021 ) »P (255.4) . V(255) ,U(255»4) 
DIMENSION  MAME<5) 

EOU I VV*LPIICE  tHtS)  ,P(  1 . 1> > • 

DATA  HftlG*<5‘»:*  * . * . D'.'flT'/ 

• DATA  X.V.W-1530*0.0/ 

DATA  lTEtF,ISTATE.ITM2/'3*l/ 

DATA  !0.  IPtl.  I PH  1 . IPN2.  ITAP/4. 255. 254,510.  *27*/ 

DATA  IPN21 » IPH4.  ICSDR/51 1.  1020.  *177570/ 

WRITER.  100) 

READ (6, 200)  MIAtEC  I) » I«l»3) 

WRITE<:6.300)  • 

REAP <6, 400) IF IL 
END  FILE  6 
IFIL-IFIL-1 
Y(  IPH) -FLOAT ( ISTATE) 

DO  15  1-2. IPH 
ISTATE«IFSR(  ISTATE.  ITAP) 

IF  (ISTATE.AI1P.  1)  5.5.10 
5 I TEMP*- 1 

10  Y(  IPII+I-1 ) "FLOAT  < I TEMP) 

15  I TEMP" 1 

DO  20  1-1.IPN1 

20  Y( I) -Y( I+IPN)  . 

CALL  LAGIMTC 10. IPM.H) 

CALL  SETFIL(3. '5T0R1 .DAT' , 1ERR. *DK* .0) 

DEFINE  FILE  3< 16. IPN2.U. ITM2) 

CALL  SETF IL ( 1 . ' TESTD . DAT* . 1ERR. ' DK' ,0) 

CALL  SETF IL<2, NAME. IERR, *MT* .0) 

DO  21  1-1.16 

21  IJRITE(3'1 ) (U( J.  1) , J"l,  IPH) 

DO  23  K-1.2 
READdXZ<J).J-l.lPN2) 

23  READ ( 1 ) (2 ( J)-»  J"  IPN21.  IPM4) 

DO  40  IK-1. IF IL. 2 

IK2-IK-2 

DO  35  IL-1,2 

IJ-IL+IK2 

J«IIQD(1J.4). 

I JJ= J*IPN 
IJK-J+1 
I JK5-J+5 
IJM-4*J 

DO  25  I -I.  IPtl  . 

rj js-ijj+i 

W(  I . 2) -ZC I JJSXP (I,2i 
UK  1.3)  -Zd.T.tS)  TP  (1,3) 

25  un.4)^(t;..i5r.ipn,4) 

I'O  26  1-1.4 

m-umi 

26  WHITE'S*  lJIXWCH,  D.Jl’l.lPH) 

IJI'IMJ 

•CALL  STASH < I Jl'N.  ICSDR) 

DO  30  K*  1 .4 


Figure  A.l  RECDll 


A-2 


LK»t10D<<lJK5-K),4)+i 

m-5-i< 

JJ2«4*<KK-1)+KK1 

REnD(3’ JJ2> (U(JJ1. D.JJ1-1.1PN) 

COLL  LnOPCIPH.X.Y.U) 

30  CONTINUE 

URITE^X’ 

if  go  ro  33 

to  32  1*1, 1FN 

32  V(I)»;;(I)St2tV(!)  . 

33  lF(UIL-l)*IIO.GE.IFIL)  GO  TO  40 
T>0  35  I-l.IPH 

35  X<I)«0.O 

IP»<I10PC1K.4)-1)*IPN 

RE0D(1)(Z(J+IP).J-1,IPH2) 

40  COtITIIIUE 

END  FILE  2 

END  FILE  1 

ENP  FILE  3 

WRITE (5.500)7 

100  FORMAT ( ’ ENTER  NAME  OF  OUTPUT  DATA  FILE’/) 

200  F0RMATC3A2) • 

300  FORMAT ( ' ENTER  All  EVEN  HO.  OF  RECORDS  TO  BE  PROCESSED'/) 

400  FORMAT ( 15) 

500  F0RIIATaX,rE10.3) 

OOO  F0RIWT(1X.5I6) 

END 


Figure  A.l  (Continued) 


A-3 


; 

; 


PEEK: 


STASH: 


LSH: 


THIS  PROGRAM  ALLOWS  FORTRAN  PROGRAMS  TO  READ  AMD  WRITE 
ABSOLUTE  LOCATIONS  III  MEMORY 

THUS  IT  IS  DANGEROUS  IF  IT  GETS  INTO  THE  WRONG  HANDS 


.TITLE  PEEK 
.GLOBL  PEEK- STASH 
.I1C  ALL  .PARAM 
.PARAH 

MOV  02<R5).R0 
MOV  <RO),R0 
RTS  R5 


MOV  02(R5) -R0 
MOV  04(R5).R1 
MOV  R6. (R1 ) 
RTS  R5 
.END 


Figure  A. 2 Function  STASH 


.TITLE  LSH 

FUNCTION  I.SH — LOGICAL  SHIFT 
CALL  IS: 

I WORD "LSH ( I WORD.  ICMT) 

I WORD:  WORD  WHOSE  BITS  ARE  TO  BE  SHIFTED 
ICIIT:  SHIFT  I WORD  BY  I CUT  BITS —POSITIVE  ICNT 
INDICATES  A LEFT  SHIFT 


.Gt.OBL  LSH 
TOT  (RS)+ 

MOV  n(P5)+.R0 

mov  cm:r5)+,ri 

UFO  PTH 

URL  If. 

rop  no 

me  pi 

ASM  P I . RO 

PTS  R5 

.END 


.•GET  I WORD 
8 GET  ICNT 

:NO  SHIFT  IF  ICIIT*U 

illO  P ROB  LEI  IS  ITU?  I NT  ' .11  IF  f 

;SH  I FT  P l SI  IT  MIKE  IT)  CI.FiTP  MSB 

:...&  COUNT  THE  SHIFT 

: SHIFT 

; RETURN 


Figure  A. 3 Function  LSH 


.TITLE  LOOP 
.GLOOL  LOOP 

ACO-X0 
l -;<0 
J*Xl 
XnD-^2 
YAD-X3 
UA  [»■?.' 4 
R5-X5 
SP-X6 

I PHI  Is  0 
YADII:  0 
UAWH  0 


LOOP: 

CMF 

CP5) . CR5)+ 

INCREMENT  R5  POINTER 

MOV 

0 (R5) +«  I 

LOAD  VALUE  OF  I PH  INTO  REG  I 

110V 

I.J 

INTO  REG  J 

MOV 

I.  IPIIM 

All!)  INTO  MEMORY  LOC  IPIIM 

MOW 

fP'51+.XftP 

LOAD  ADPR  OF  APPAY  M IHTO  PEG 

XAD 

1 MV 

(fiSJ+ViTID 

LOAD  HDDR  OF  Al  l AY  III  III  l:i-.G 

vnii 

SLID 

*4. YAP 

SUBTRACT  4 TO  POINT  TO  PREVIOUS  IJORD 

mv 

Yf'tD.  VADM 

MOV 

CRS)+,IJAD 

MOV 

UhDAIADM 

MOV 

R3. -<SP) 

Lis 

MOV 

I.R5 

ASL 

■ R5 

ASL 

R5 

ADD 

R5.YAD 

LUF 

( XAD ) , ACO 

L2: 

MOV 

(Yi'iD)  +.R5 

BUI 

SUBTP 

ADIiF 

(LlnrO+.AC0 

EtR 

.+4 

SUBTP 

sSUBF 

mjaio+.aco 

TST 

Ci'hD)  + 

DEC 

J 

RUE 

12 

, 

STR 

ACO  . C<AD)  + 

1 10V 

IPIIM.  J 

1 Kiv 

I'ADII..  WilD 

■MOV 

Yl'iDI  I,  YAD 

DEC 

I. 

01  IE 

LI 

1 10V 

(SP)4-.R5 

RT& 

R5 

.El  ID 

Figure  A. 4 Subroutine  LOOP 


i — MTIN. MAC  — 

; PURPOSE  IS  TO  REAP  III  FROM  MAGTAPE 
i CALL  BY  COLL  HrilKlHIlT.  lWORPS.  IiTPPAY.  ISTAT) 

WHERE  IJHIT  « O OR  1.  I WORDS  13  NUIIBER  OF  WORDS  TO 
TRANSFER i 10O0-).  1 nf. ‘PAY  IS  THE  ARRAY'  THE 
IHFO  IS  TO  GO  INTO.  PHD  ISTAT  IS  A STATUS 
WORD  PlilUPIir.il  WITH  FOLLOWING  MEANINGS: 

ISTAT  - 0 •>  O.K. 

» 1 ->  EOF  FOIJIID 
; - 2 -■>  UNIT  SELECTED  IN  ERROR 

THE  MAC.  TAPE  STATUS  WORD  CUTS)  13  RETURNED 
IF  THE  EPROR  BITS  (T- 155  OF  (ITS  ARE  SET 
*"M*1  WORDS  MUST  BE  AH  INTEGER  VARIABLE 
AND  HOT  A CONSTANT  AS  INFORMATION  IS  RETURNED 
IN  THIS  VARIABLE  AS  TO  HOW  MANY  WORDS  ARE 
TRANSFERRED. 

.TITLE  ITT' I H 
.GLOHL  I1TIN 
. I CALL  .PARAM 
. PAPAM 
INS- 172020 
IITC"  172522 
I ITBRC- 172524 
I ITCMA -172526 


i PO-UNIT  * 

;SET  BYTE  RECORD  COUNT 
i R2 -ADDRESS  OF  ARRAY 


MTIN i MOV 

RO.-(SP) 

MOV 

Rl.-(SP) 

MOV 

R2.-(SP) 

TST 

(R5)  + 

MOV 

|J(R5)  + .P0 

MOV 

|J(,RS)  + .R1 

MOV 

(R5)+.P2 

ASL 

R 1 

■ HEG 

PI 

MOV 

R 1 , l?*HTBRC 

BIT 

ttirrrrc.po 

BEO 

COHT 

CLP 

0-4 (P5) 

MOV 

«2.0(R5)4- 

BR  OUT 

; TEST  UNIT  <1  OR  0) 


COIIT:  MOV  P2.GMN1TCMA 

SWAB  PO  i 

ADD  ttG0003.R0 
TTs  MOV  ttllTC.Rl  ; 

TSTB  ( P n ; 

PPL  .-2 

1,10V  RO.  (PI  i : 

. PAGE 

MOV  III  I IS . hi 
Wli  DIT  112000. ' PCO 

DEO  1,12 

CLP  (PI)  : 

uii!  W3 

W2:  TSTFj  'PH 

FJPL  l.l  I 


; NO  WORDS  TRANSFERRED 
; ERROR.  RETURN 


SET  COMMANDS 

R 1 - ADDRESS  OF  MTC 
TEST  CU  READY 

SThRT  the  operation 

i EO" AD PRESS  OF  MS 
: HALT  IF  EOT 

STOP  Ti-lPE 

TEST  CU  READY 


Figure  A. 5 


Subroutine  MTIN 


Piirnsiou  x<,’ooo).v(*iio) . iop<  inon:> 

DATA  11.112.  IS. IP.  (OFF  lOOO.2nnn.n.51O.0>' 
Ci'iLL  ->ni'l  II. ' 1.  * 1ESTL'.  I'll  I ■ , 1ELR.  * bi:  * .0) 
t.rj7F<fi.  f5ft> 
h irsiLL 

' IFCLL.EO.OGfl  TO  00 
Ki  ?'.»  H.LL 
DO  TO  J«1 .3000 
I WORDS* 1000 

CALI.  MTIIIfO.  II.I0RPS.1AP.  ISTAT) 
IFdSTAT.EP.  HGO  TO  75 
WRITE <6.  175)  J 
WRITE (b. 125) 

PEAIK6. 1 75)  LL 
I WORDS -1000 
IF 'LL.EO.OIGO  TO  3 
DO  I 1-1  - H. 

C ALL  HTimn.  IUOPPS.  1AR.  ISTAT) 
IFdSTAl.KO.  1*1.50  TO  200 
IFf  J'.mPDS.MF.  innO)GO  TO  225 
DO  2 >1.2 
j-n-m 

CALL  tiriHLO.  1 WORDS. 1AR, ISTAT) 

IF(  ISTAT. EO.  1>G0  TO  200 
IF (1 WORDS. HE. 1000)  GO  TO  225 
DO  2 1*1 .H 

:■  :c  1 +J 1 1 *FLOAT(LSH<  I HR  ( I ) . -4)  -*4000) 
WRITE (6/275) 

PEAKS,  ITS*  IF ri. 
lFdFIL.GT, lOOO) IF IL- 1030 
DO  40  J*l,  IF  II. 

I F < (112-  I OFF)  .LT.  IP)  GO  TO  15 
DO  10  1-1,  IP 

Yf  !)»>:(  I + 10FF) 

I OFF"  I (IFF  IIP 

iFfris.no. n .of;. uoff.lt. in)  go  to  30 
COLL  NrillfO.  IUOPPS.  IAR,  ISTAT) 
IFdSTAT.EO.il GO  TO  200 
I F d WORDS . 1 IE . 1 000 ) GO  TO  225 
DO  12  1*1,11 

):<  n»FLOAT(LSH(  lARd),-4>-"4000) 

IS-1 

00  TO  30 
I J -M.:~  toFr 
DO  20  1*1,  I J 
Yd)-Xf  I + IOI'FJ 
113"  1 

n" ' I i .r.r.m  ii>2 
DO  .•'?  ! r I = 1 . Ii; 

CHI  L Ml  III  CO . I WORDS.  I IIP,  ISTOT.I 

ii;'.  mm  .r:n.  noo  to  200 
) F 1 1 i .'iJI  'DS  . I IF . 1 000 ) GO  10  225 
ill  'l  l 1 j ••  1 :*  hi 
IF  1 113. EO.  1'  H>H 


Figure  A. 6 Subroutine  FORMlA 


TCI  22  I- Ml 

22  >:<  l+Nl ) -FLOAT (LSH(  IAR(  I)  .-4) -’4000) 

IS-0 

IJ1-IJ+1 
I'O  .?*5  I -Ml.  IP 
25  Y(  1)  »>'(  J-IJ) 

iorF»iP-u 
30  l.:p  1 1 E k 1 ) V 

WRITE'S.  IOO)J 

IP  ( ( 10FF  . LT. II)  .OR.  < IS.EO.  1) ) GO  TO  40 
COLL  (ITlHfO.  I LORI'S.  IrtR.  ISTAT) 

IF ( ISTAT. EO. 1)60  TO  200 
IF  C I WOW'S.  I IE.  1000)60  TO  225 
I'O  35'  I*I.N 

35  X(I»-FLfjnT(I.SH(IAR<I).-4)-»40O0) 

1**1 

10  CPI  IT  I HUE 

l no  fopiiatukiois) 

175  FORlirm-  FIITEP  IIP.  OF  RECORPS  TP  BE  SKIPPED'/) 

I'.ij  FORIIM  I < ‘ I'll  I CP  IIP.  or  P ILL'S  TO  DC.  SKIPPED’/) 

IPS  rpPIVlTCHv' 

200  EMI*  FILE  I 

STOP 

22*5  WRiTi?i.e.p*j:n)  rvmT 

250  rCiRIIOTr*  PI  i'll1  ERROR’/’  MT  5T0TU3  LIOPT):  ’.06/) 

2r’5  FORI l:;iT ( ’ U I I EP  MO.  OF  RECORDS  TO  BE  PROCESSED’/) 

CHI'  F ILE  I 
Elll' 


Figure  A. 6 (Continued) 
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.TITLE 
.HL 1ST 


SDEVTB 
TTM.BEX, TOC 


;THIS  IS  A STRIPPED  DEVICE  TABLE  FOR 
j NRL  • 0 

NPLI1T  - 0 

JNRLMT  - 0 


USE  WITH  PLAY1 1 ONLY. 

••IF  ' HPL ' IS  DEFINED  USE  NPL  DEVICES 
; IF  URLMT  IS  DEFINED  USE  (IT  AT  NRL 
MF  HOT  DEFINED  USE'W'l:  AT  NPL 


*DEVTB: : . WORD 

3 

. 

.LORD 

0.DEV1.PEV2.DEV3 

DEVI:  . LORD 

0 

: DDB  POINTER 

. RAD50 

-■DKV 

.•UNIT  l DEVICE 

.BYTE 

0 

SHOW  CURRENTLY  OPEN 

.BYTE 

0 

.•UNIT  1 DEVICE  UNIT  NUMBER 

.RAD50 

/PLAYD/ 

: UN I T 1 FILE  HATE 

.RAD50 

/PAT/ 

;UN1T  1 EXTENSION 

. UOPD 

233. 0»0. 0.0.  ,0.0. 0,0,0 

:OTHER  USEFUL  INFORMATION 

DEV2 : ■ . LORD 

0 

. PAD5Q 

,111/ 

Jim  IT  2 DEVICE 

.BYTE 

0 

' .BYTE 

0 

JUNIT  2 DEVICE  UNIT  HUMBER 

.RAD50 

/PECDT/ 

JUIIIT  2 FILE  MAI  £ 

. RAD50 

/DAT/ 

.’UNIT  2 FILE  EXTENSION 

.LORD 

233.0.0.0.0.0,0.0,0.0 

, 

DEV3:  .LORD 

0 

. IF  NDF  NRL 

; IF  'NRL'  IS  NOT  DEFINED  USE  THIS  FORM 

. RADSO 

/DK/ 

. BYTE 

0 

.BYTE 

0 

;UNIT  3 DEVICE  UNIT  HUMBER 

. IFF 

.IF  DF  NR1.MT 

; IF  ’URL'  AND 

‘IIRLIIT'  ARE  DEFINED  USE 

THIS  FORM 

. PADS0 

/I  IT/ 

. BYTE 

0.1 

. IFF 

: IF  'NPL.'  IS  DEC  HIED  AMD  'IIRLIIT'  IS  NOT  DEFINED  USE  THIS  FORM 

.PAPSO 

- J'l  •' 

. Hi  IE 

1 

. ENDC 
.ENDC 

.RAD50 

TROOP/ 

JUIIIT  3 FILE  tlfll IF 

.T’tiDLiO 

••  DAT  ' 

JUIIIT  3 EXTENSION 

. DARI) 

<i’ 5 ^ * 0 0 0 * 0*O*0*0*0*0*0 

. END 


Figure  A. 7 Device  Table  MINDEV 
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.TITLE  ERRF . . 

; MINI HAL  ERROR  HANDIER  FOP  FORTRAN  OTS  V22 

: OUTPUT  MESSAGE  HA?  IDENTICAL  INTERPRETATION  TO  FORTRAN  OTS  GENERATED 
; F030  ABCKYZ  MESSAGE.  ABC-EPROR  CLASS.  J^YZ -ERROR  NUMBER 

JHOTE  THAT  EVEN  THE  MOST  TRIVIAL  ERROR  DILL  BE  FATAL 


ERROO0:  :110V 

O2(R5).R0 

DR 

-TERRA 

TERR  s : MOV 

2CR5).P0 

JERRA : : 

*ERRB::  SWAB 

R0 

1 10  VB 

R0.RI 

sue 

PI.R0 

A SI. 

PO 

Hl'D 

r i .pci 

1 10V 

F'O.  - fSP) 

1 10  V 

#M30.-(SP) 

IOT 

EHT 

60 

MTRRF  i : 
■£ERRC.S ! 


.END 


Figure  A. 8 Subroutine  ERRF 
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CONIOIVOVPxYUZe^.ZC  lO20).W(12e2),H(1031)  ,.X<  1020) 

COUNOH.'PAVAL.'II.  ICSDR,  IQ.  ipn.  ipmiq,  ipimqi,  iset. siphon 
DIMENSION  P(  1020) 

EPUIVALEHCE(H(2>,P<1>> 

IPUKl  “10  f I PI  I 
IPIIIOl-lPHIO+IPH 

SIPIID4“1.0.'i4.O+FlOAT(  1PN+1)#12O0.) 

IREC-3 
uni-i 

call  Lonriri.AG' . 1.  IER) 

FILE1:-.  APE  ASSIGNED  WITHIN  THE  IIIIIIMUM  DEVICE  TOOLE  PPOGRAM-MINPEV. 
IF  THAT  FPCIGPAM  IS  REMOVED.  THE  FOLLOIJIHG  THREE  CALLS  TO  SETFIL 
ARE  REQUIRED . 

COLL  RETF 11(1. 'PLAYD.DAT' . IERP. ' I«K' , 0) 

CALL  SClFIL(2,'PECI>T.PAr>  1ERR.MIT' .0) 

CALL  SETF IL (3. ' PROCD.DAT' . IERR. ' DK' • 0) 

CALL  LIAIT 

CALL  LAG IMT ( 10. 1PM. H) 

CALL  LOAD (‘SUB*. 0. 1ER) 

6 CALL  SUni(IPEC) 

IF1.  ISET. HE. COCCI  TO  GO 
IJK'OI  Jl'N+i 

CALL  STASH ( l.M(H. ICSDR) 

CALL  LOADCFFTOVR'.O.  IER) 

CALL  FFTOVR 

CALL  LOADCSIJB',0.  IER) 

IF  CUKfl.  HE..  Q)GO  TO  30 
II  — I 

CALL  SUB2UPN) 

II -3 
1PEC-] 

GO  TO  6 

30  CALL  SUB2(0) 

GO  TO  6 

SO  Fill)  FILE  3 

clip  F ILE  3 
FMP  FILE  1 
1 .11  1 1«  I. II  1 1 1 I'M:  T 
C ALL  AT hSI I ■'  1 ll-'H.  ICSI'P) 
till' 

BLOCT  DATA 

torn  toil  OVR.  V 1202) . E C 1020) . IJ(  1202) , H(  1021)  1020) 

COHIIOH  I'AVAL.  f I.  ICSDR.  Ip.  IPII.  IPHIQ.  IPHIO  1 . ISF.T.SlPWM 
DAI A V.2  230210. O' 

DA  I' A I I,  ICSDR.  IQ.  IPII.  ISET.  DO."  ITTSPO. 4,255.0.' 

PHD 


Figure  A. 9 PLAYll 
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SUBROUTINE  SUBIC IREC) 

COItlON.'OVR/'YC  1 282).  Z ( 1020),  UC12a2),H  002  n.MC  1020) 
COIUiN.'DNVhU'II.  ICSDP.  10.  IPN.  IPII10.  IPNIQl.  IBET.SIPND4 
1-0 

m-c  i+i  duipii 

PE8DC 1 .END’S) CXCJ+II1),  J«I, IPH) 

I - I + I 

IFCI.LT.  IREOGO  TO  1 
IFC  IREC.GT.  DGO  TO  4 
1*0  2 1-1 . IPH 

vu+iniitn-xn+im 

do  3 t-1.7 

yi  i +ipiuri i > "0 .o 

PEmI'<2.EMI-SUJ 

RETURN 

ISET  " "‘10000 

RETURN 

I SET- "100000 

RETURN 

END 


Figure  A. 10  Subroutine  SUB1 


SUBROUTINE  2t.IRITECZ.IPN) 

PEmL  z c jpi-o 

UP  1 TECS) z 

RETURN 

END 


Figure  A. 11  Subroutine  ZWRITE 


.TITLE  SU02. 

.ML  1ST  TTM.BEX 
, PSADL  GBL 
, GLOBL  ZURITE 
AC0  • y.Q 

AC1  « t- 1 

j MACRO  REPLACEMENT  FOR  SUBROUTINE  SUB2 


•CSECT 

Cl  VP 

Y: 

.BLKU 

1282. *2 

Zs 

.BLKU 

1020. *2 

U: 

.ELI:  U 

1282.12 

H: 

. BLKU 

1021. *2 

Xt 

.BLKU 

1020. *2 

P 

■ 

H+4 

.CSECT 

DAVAL 

III 

.BLKU 

1 

ICSDR: 

. BLKU 

I 

IQ: 

.BLKU 

1 

I PH: 

.BLKU 

1 

1PM 10: 

.BLKU 

I 

I PH  101: 

.BLKU 

1 

ISET: 

.BLKU 

1 

S1PND4: 

.BLKU  ' 

2 

; FORTRAN  COMMON  REGION 


; ANOTHER  COMMON  REGION 


.CSECT 

. SUB2. 

III: 

.UORP 

0 

SUB2: : 

MOV 

II»R0 

MOV 

SETF 

R5» -<SP) 

INC 

P0 

CMP 

PO.  10 

BMP 

1$ 

CLP 

RO 

11: 

MOV 

PO.  11 

IIUL 

IPN.R0 

MOV 

PI.  11 1 

BGT 

34$ 

TST 

02(R5> 

BBQ 

31$ 

29$: 

MOV 

IPII.RO 

ASL 

RO 

A3L 

PO 

riT-ii 

ikV.  PO 

1 10 '/  . 

1 10  V ' • 

IFIII0.P2 

30$: 

1 10V 

i.l'IIL  >.R0)  + 

1 IOV 

iRl >+. (RQ)  + 

son 

P2 » 30$ 

BP 

55$ 

Figure  A. 12  Subroutine  SUB 2 
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77*1 

LDF 

(RD+.AC0 

MULF 

<P2)*,AC0 

IIOV 

(R4)+. <R3)  + 

ADDF 

(RO).ACO 

MOV 

<R4)+, <R3>+ 

STF 

ACO» <R0)  + 

SOB 

R5.37* 

384: 

MOV  ' 

11I,R0 

11-r.L. 

ro 

ASL 

F!0 

ADD 

*Z«R0 

MOV 

R0.43* 

MOV 

R0».-(SP) 

MOV 

IPN.R5 

MOV 

R5.-CSP) 

LDF 

SIPND4.AC1 

394: 

LDF 

(R0J.AC0 

MULF 

AC1.AC0 

STF 

ACO,(R0)+ 

SOB 

R5.39* 

JSR 

R5f ZU?ITE 

BA 

44* 

434 1 

.WORD 

0. 1 PM 

4441 

MOV 

(SPJ+.R4 

1 10  V 

(SP)+.R1 

,45*! 

CLRF 

(R1H 

SOB 

94. 45* 

554: 

MOV 

(SP)+.R5 

RTS 

R5 

.END 


Figure  A. 12  (Continued) 
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SUBROUTINE  FFTOVR 

COIHON/OVR/Y(  1282)  ,2C  1020)  ..u<  1282) >H(  1021) ,XC 1020) 
COMWH/DAVRl/'II,  IC?t>R.  10,  IPN,  IPNIQ,  IPMIQl,  ISET.SIPND4 
COUPLE):  YC<6-11).l.rC641) 

EOlJ  1 VhLHIICE  CY<  1 ) . VC H ) > . OK  1 ) , UC  ( 1 ) ) 

COLL  LOAtirFFT' ,0.  IER) 

Cf’ll.L  rnrvt,  I ) . V<2), 640. 640,640. 2) 

COLL  LOOt't'REOLTR'.O,  IER) 

COLL  RGflLTRCYf 1 ) ,Y(2) >640,2) 

£•0  2 1-1.641 
YCU)«yC<I>HJC(l) 

COLL  REOLTR<Y(l).Y<2).640.-2) 

•COLL  L001K  *FFT'  .0.  IER) 

COLL  PFT(V<  I ) ,Y<!2)  .640,640,640,  -2) 

RETURN 

END 


Figure  A. 13  Subroutine  FFOVR 
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c-riLf!"  muvwTT  , csio.r'.ioi  04-jm-74 

uidpoijuni-  fftca.b.htot.n.hspnh.  rsm 

I' II  If- MS  I Oi  l rnO.B(l) 

PirCHSlOM  HFi1C<0).NP<9) 

I> II ICIIS ION  riT(5).BT(5) 

EOtm'filEHCE  U.1I) 
rnxF-s 

iri:\P-9  . 

l»0 

JF-0 

K*  I 

M"8 

KT-3 

HFflC<t)-2 
HFAC(2)»2 
NFAC<3)-2 
NFAC<4)-2 
IIFAC^-S 
HFrtC(S) »2 
NFAC<?) “2 

l{FAC<e>-2 

IIIC-ISII 

ftAD-B.OnVTANU.B) 

S72-RA1V5.0 

Cr2«CO&tS?2> 

$?2«SIH<S?2> 

Sl20»SORT<O.fi5) 

IF  <ISM.GE.0:i  GO  TO  10 

Sr2»~SF2 

S12O--S120 

PAS-RAD 

INC  — II  1C 

0 ht- inc  h hot 

KS-IHCHISPAH 
k'SPAN«|:S 
UN-NT -INC 
JC-KS-N 

RADF-PAMFI.  0AT(JC)>i;.5 

ao  OD»RnDF.-Monra:spnH) 

Ct**2 . t'  i"S  I N ( SIO  Hcftg  ■ 

CD-SlN'SD+SD) 

M>1 
I ■ I > I 

IF  nir-NC  I)  .NE.2)  GO  TO  600 

KSPnn-i  nPnn^2 

Kl-ly:iFTillt-2 

io  i-'2*ti'+i.:,phn 

[il  >1.1 1 1 Z’  l 
Nrk2i  ■ h I W I -I  II: 

B 1 1:2  I "fji.l.i:  I f.1K 
fl  rIT>  ••|‘l 1 l;l"  > INI; 

Bi  i i.''r.iii  r 1 4BI-: 


Figure  A. 14  Subroutine  FFT 


IT  • li.Mi  Sfnll 

jr  *1  ll.LC.llll>  C.0  TO  210 

l»i:K-MH 

IF  (I  I .LE.JC)  (.0  TO  210 
IF  m..GT.KSPAII)  GO  TO  800 
220  Cl-I.U-CB 
S1«SI> 

230  K2’KK+KSPAH 

AK*A(KI)-A(K2) 

BK*Bd  K)-B(K2> 
A(KK)»A<KK)+A(K2) 
OO'.IO-B(KK)  fB(K2) 
A(K2'«Cl*rtK-3tl15K 
B <K2) »S 1 iTiK+C  l-'  BK 
KK-K2+KSPAII 

IF  (KI'.LT.IIT)  GO  TO  230 

K2-KK-IIT 

Cl— Cl 

KK«K1-K2 

IF  (KK.GT.K2)  GO  TO  230 
AK«CI-<CIdCl+S!»»Sl> 

Sl"(SI»  »U-CI>*SI)+S1 

Cl-AK 

KK-l'K+JC 

IF  (KK.LT.1'2)  GO  TO  230 

KI-KI+IIIC+1HC 

KK  ■ <K  1 -KSF'AM)  /2+JC 

IF  d'K.LE.  JC+JO  GO  TO  220 

GO  TO  100  • 

510  C2-C72+*2-$72*'i2 

S2-2.0tC72i572 
520  Kl-KI+KSPAH 
K2-K1+KSPAH 
K3-K2+KSPAH 
K4-K3+KSPAH 
AKP»A'K1>+A(K4) 

AKt1*A<K  1 )-Ar(K4) 
BkP»B'Kl)+B(K4) 
BI<H*B(K1>-B(K4) 
AJP«A(K2)M)(K3) 
AJM*Ad:2>-Ad'3> 
B..tP*P(l'2MB(K3> 

BJN'Bd  2>-B(K3) 

AA»A(I'!K) 

BB-B(KK) 

A (IT.)  “AAtAKP+AJP 
B(l.l ) -BB+BKP4BJF 
Al : =h'-  P i C 72+h  IP:I  C2+AA 
BK  *BI'.P  ^ C72+0  JPIC2+BB 
A J “O'- 1 1 ' 372+A  M I > S2 
B.T=BIHi572»BJIIi«2 
h«ki:<  -f'  ,-BJ 
Ad.  4 1 -AK+B J 
BMP  -"PP  »0J 


Figure  A. 14  (Continued) 


fii:*oi;P*t24njptcr2+ftn  . 
Bi:«nr.p  iC2<njnc72iBG 

BJ«ei.tHS2-BJN»S72 

BtW-'AK-BJ 

H(K3) 

E!iK3’i  =BK-flJ 

KK"K4+KSPrtll 

IF  (KI..LT.IIM)  GO  TO  520 
KK-KI.'HIM 

IF  (KK.LE.KSPrtH)  GO  TO  520 
GO  TO  TOO 
600  IMIFGCn) 

KSPHtl-KSPOH 
KSPAH*rSPAM/K 
IF  (K.EQ.5>  GO  TO  510 
700  IF  (I.EQ.ID  GO  TO  800 
KKrJC+l 

710  C2* 1 .O-Ct- 

SI  •St- 
720  Cl  *C2 
S2-S1 

KK-l'.K+KSPflH 
730  AK'fiCMO 

AvKIO-C2-'»AK-S2*0<KK> 
B<KIO-S2*M'+C2*B<KIO 
KK-KI.  H.SFIIH 

IF  (KK.LE.HT)  GO  TO  730 
AOS1+S2 
S2-SI «C2+C1*S2 
C2-CI i C2-0K 
KI:-i:i>HT+KSrOII 
IF  (KK.LE.KSPMH)  GO  TO  730 
C2*C  I -•’CMC  H-$I»*S1> 
Sl-SIMSWCI-CD+Sl) 
ki:-ki;-i;sphh+JC 
IF  (kK.LE.KSPHH)  GO  TO  7 20 
i:i:*i:t;-i'SPOHfjc+iHC 
IF  (KK.LC. JC+JC.'  GO  TO  710 
GO  TO  1 00 
300  IIP  ( 1 - =KS 

IF  (I.T  .-E'1.0)  GO  TO  890 

K=kT+KT+l 

IF  (H.LT.IO  l>K-l. 

J»1 

HP  - 1:  4)  • 

3io  iiPt.i  h > -up-  j.'  iirnc(J) 
i ip 1 1 ■ -up -I  t-r-  mi  nc<;j) 

J“.M! 

► «)■:- 1 

IF  (J.LT.K)  GO  TO  310 
K3-MP0+1  > 


I l*i ill- HP i .?) 

I I • .li.  i- 1 
I '2  * I SPirt If  1 
J-  1 

AKWHI.K) 
aim:j  i:2> 

t’i(N!2)  "i  iK 
DI>OCKIO 
B<KI.->»B<K2) 

EUK2>»BK 

kt-ick+iiic 

k‘2*K«.PAII+K2 

If-  (K2.LT.ICS>  GO  TO  820 
K2»K2-Mri.,l) 

J-JM 

l<2,,IIP<JH)’ll<2 

IF  <K2.GT.NP<J>>  GO  TO  030 

J-l 

IF  (KK.LT.K2)  GO  TO  820 

kk-mc+ihc 

K2-l'5PAm'K2 

IF  (l:2.LT.KS>  GO  TO  040 
IF  (KK.LT.IC5>  GO  TO  830 
JOK3 

IF  C. Hi-: T-H  . GE . M)  RETURN 
KSP,,H“HP(KT'H> 

J-IMCT 

HFAC(J+»)«I 

IIFAC<J)»IIFAC(J)KIIFFiC<J+l) 

J-J-l 

IF  U.HE.KT)  GO  TO  900 

KT-KT+I 

HH»HFfiC(KT>-l 

IF  (HH.GT.HAXP)  GO  TO  998 

JJ-0 

J«0 

GO  TO  906 
JJ-JJ-K2 
IC2-KIC 
K-r+i 

KK-IIFACOO 
j.i»fcr  Ki..i 

IF  ( JJ.GF. .I:'2>  GO  TO  902 

UPC Ji - JJ 

K2*HFflC<ICT) 

k'-KT+l 

KIC-  HI  i'iC  CIO 

J-J-M 

IF  (J.LE.IIH>  GO  TO  904 
J»0 

GO  TO  914 
l»kt. 

I K '-I  ip  v I : ,i 

HP* l:  I —IK 


Figure  A. 14  (Continued) 


IF-  ikl..llfc\J)  GO  TO  910 

1.3-1 1: 

J-JH 
l:l'“IIP(  J) 

IF  (i:l'.I.T.O»  no  TO  914 
IF  •IKU'.IIE. J>  GO  TO  910 
HP  < J :•  —J 

IF  (J. ME. Nil)  GO  TO  914 
MrtXF.  Hh’FiCIMC 
GO  TU  950 
J'J-I 

IF  1 1 IP  ( J ) .l.T.0)  GO  TO  924 

JJ-JC 

F5PWKJJ 

IF  MJ.GT.IMXF)  KSPAIWWXF 

JJ-JJ-KSPRII 

K»Hr>(  J) 

KK-JC  I i:-»  I I+JJ 
K 1 -KK+KSPAH 
k2»Q 
K2-U2+1 
rtT(L2)»ft<KI)  • 

BT(i:2)-B0<l) 

Ki-i.i-mc 

IF  (KI.IIE.KIOGO  TO  928 

li-KK+KSFAH 

i:2-l:  1-JC*<K+HP<K>> 

I>.-NP<K> 

ACK  l W1CK2) 

BiKl)-B(K2> 

K 1 -K 1 - I HC 
K2-K2-INC 

IF  (K1.IIE.KK)  GO  TO  936 

1 i:-i:2 

IF  (K.IIE.J)  GO  TO  932 

k 1 “KTMISPON 

1.2-0 

I',.-?' IM+l 

H (Kl)»ftT  <T2'' 

tJ'.K i ■»  -I3TCK2> 

i 1-1.1-  inr 

IF  'T.1.HF..KIO  GO  TO  940 

IF  - JJ.MF.O)  GO  TO  926 

IF  J.IIE.  1)  GO  TO  924 

J-K3+1 

t IT -MT-kOF'HN 

II-MT-1NC+1 

IF (NT.GE.0)  GO  TO  924 

RETURN 

ISH-0 

CALL  ERR0O0(265> 

Fill' 
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SUBROUTINE  REfiLTR(fl,B<N» ISN) 
I) IICNS ION  A(  1 ) ,B( 1 ) 

PEAL  III 
IIIOIAGS<ISN) 

ni.-ii  miiic+2 

NIWIIV2 

SD*,?.0+rtTi1ll(  1 .0)/FLOAT(H) 

CI>'2.04SimSD)**2 

pr'-siiusc+SDi 

SH«0.0 

IF  CISN.LT.Q)  GO  TO  30 
CM-1.0 

0(Mk-D«B(l) 
t'O  20  J-I.HH,  INC 
K-NK-J 

Cirt  “D  < J ) -t-B  1 . 

PB»EUJ)~B(IO 

PE-CIMOrt+SH'(flB 

IH'SN*BM-CH*rtB 

D(IO-IllrBB 

EKJ)«IN+BG  1 

J) -ftA+RE 

AN-CIWCt>+CM+5D*SH> 

SI  1“  ‘.SI'iC  H-CD*SN>  +SN 

CH»AA 

RETURN 

CN— 1.0 

sr>— sd 

GO  TO  10 
EMI) 


Figure  A. 15  Subroutine  REALTk 


SUBROUTINE  SUB2<IOFF) 

COMMOH/OVP/Yt 1202) ,2( 1020) .W< 1202) .H( 1821) .XC 1020) 

COI I WH-  DAVAL  '11.  ICSDR.  IQ . IPH.  IPN 10.  IPH 10 1 . ISET. S l PND4 
DHEH9I0H  P( 1020) 

EC!UIVALEHCE<H(2).PU)) 

I I*MQD< ( 1 1+1) . 10) 

111*1 1* IPH 

I F ( 1 1 1 ,GT.0)GO  TO  34 
IF(IOFF.EQ.O)GO  TO  31 
DO  30  1*1. 1 PM  IQ 
II2-1+1PN 

30  Y< 1 12) *X( I) 

GO  TO  53 

31  DO  32  I-1.1PHI0 

Z<  l) *Y< I+IPH)*P( 1)+Z( I) 

32  Y<  I ) *X<  I ) 

GO  TO  38 

34  II4-IPN10-III 

DO  35  1-1. 114  • 

1 13*1+111 

Z(tli)"Y<l+lPN>*P(I>4Z(!I3> 

35  Y( 1) *X( 1 13) 

I12-IPHIQ1-1I1 
DO  3?  1-1. Ill 

1 13*1+1 14 

2(1) *Y( l+II2)#P(I I3)+Z( I) 

37  YUI3)-X(I> 

30  DO  39  1-1. IPH 

1 14*1+1 1 1 

39  2 ( 1 1 4)  *2 < 1 1 4)  >kS  IPND4 

WRITE (3) (2(1+1 1 1) . 1*1. IPH) 

45  DO  50  1*1, IPH 

30  2 ( l+I 1 1) *0.0 

55  RETURN 

END 


Figure  A. 16  Subroutine  SUB2 
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.NAME  LAG 
.NAME  SUB 

i . T-XTR  PLAVll/tC-ERRF-MINDEV-CNRLIB-FTNLIB 

! .FCTR  LUG-LAG 1-FTHL  IB 

! .FCTR  bUB-SUB  l-SUB2-i'URITE-FTNL IB 

i .FCTR  FFT-FTNLIB 

i .FCTR  REALTR-FTNLIB 

: . FCTF:  FFTOVR-FTNLIB 

♦ROOT  A-<B.C. (F-(D.E))) 

.END 


Figure  A. 17  Overlay  Descriptor  PLAYll.ODL 


A-25 


j 'i  in 1. 1 -it* 


DIMENSION  ,':(255).Y<5O9).ZU0S6).M<5n),P(255.2).V:255).UC255.2),HT 
t (22 1 . HZ ( 1020 ) . ST<  46 ) . NAME  (5)  .Ffll’E  ’ 2048 ) 

PROGRAM  MUST  PE  COMPILED  IJITH  THE  Oil  SWITCH. 

EOU I VALENCE  (H  (2?  .PM.D).  (TEMP.  IAR) 

DATA  ’ .D' . 'AT'  ' 

DATA  M.V.lil.  Z.  ITFMP.  I ST  ATE.  ITH2/'2OR6*0.O.3*1/ 

DATA  I PH. IPHI.IPH2. ITAP.255, 254,510. "2“0/' 

DATA  IPII22I.  IPH222.  IPII4,  IPM421 . ICSDR,  ISTOR/53 1.532,  1020, 1041. " 1775 
8,70. "200/ 

DATA  IPII247.  IPN446/557, 1056/ 

DATA  HT/.6242G603E-2.O . . . 75582797E-2, 0 ...  1 I375S08E- 1,0.,. 1GO07 142E 
8 - 1 . 0 ...  26 1 709O6E- 1 , 0 . . . 373 1 1 5S 1 E- 1 , 0 . , . 5298709 1 E- 1 . 0 . . . 76534B23E- J 
*,0.. . 1 1669172.0. . .20568788.0. . .63442332.0./ 

IR3TAT-0 
UPITEC6, 100) 

RE AD (6. 200)  MIAME ( I ) . 1-1.3)  ' 

URITEC6.300) 

READC6.4O0) IF1L 
END  FILE  6 
IF  1L  " IF  II.  - 1 
V IPID-FLOnrUSTATfi) 

DO  11  I »2. IPN 
ISTATE-IFSR( ISTATE. ITAP) 

IF < I STATE. AND. 1)3.3. 10 
ITEMP-- 1 

Y( IPH+I-1 ) "FLOAT ( l TEMP) 

1 TEMP-1 

DO  12  I -1 . 1PM1 
V< I ) -Y( I+IPH) 

CALL  SETF IL( 1 , 'LAGR.DAT' , IERR, 'DK' ,0) 

PEAD(1)H 
END  FILE  I 

CALL  SETF IL<3. 'ST0R1.DAT' , IERR. 'DK' ,0) 

DEFINE  FILE  3(8. IPN2.U, ITM2) 

CALL  SETFILU.  'TESTD.DAT'.  IERR.  'DK',0) 

CALL  SETFILU. NAME.  tERR.'IIT*  .0) 

DO  13  1-1.8 

WRITE  (3*  I TILT)  <UC  J . I ) . J • 1 . IPN) 

READ 1 1 .END-75)  ■'Z(J).  J-22,  IPN221) 

DO  50  IK-1. IFIL.2 
IK2-IK-2 
DO  40  1L- 1,2 
U-IL+IK2 
.I-MOD(  I J . 4) 

1 .1.1  -J''!  PI  1+21 
1SI-  f— 4 

lF(.r.HF.nvr,n  to  ?3 
i irr.FTiim  rrr;i  jl.  hiring 
Du  i'.,  i-i.  in  it 

1 1 I r ll'll. I 
:: ' i Hi  '■'(  1 11-25) 

DO  16  1-1.46 
Z'P-STU) 

READ I .EHIi-75)  (ZCJQ) . J0-IPII247.  IFH446) 
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t'O  19  I "4? . IPII446 
TE I II 1 - ~ < n+Z(  1-4' 

IF < iHP.Kn.m&O  TO  19 

C SUBTRACTING  "200  EFFECTIVELY  DIVIDES  BY  2. 

IAP-IhR-ISTOP 

19  2(1-4) -Z<l-2>-TEH» 

I'O  20  1-1.  IF. 

20  ST i I • -Cl  M ini  I> 

C 'l  lCni.CUI.flTE  HILBEPT  TRANSFORM 

DO  22  I-l. IPM4 
TEMP *0.0 
DO  21  KL- 1.21.2 

2 1 TEMP -(2(1 +KL- 1) -Z ( I -KL+43) ) *HT(KL ) +TEMP 

22  HZ (I) -TEMP 

C xofORPELATF  AMD  INTERPOLATE 

23  IJIWWt»<J.2) 

IJK3-IJM+3 
IJ11-L.SH(  IJM.  1) 

24  ISET - ISET+4 

IF < 1SET.EO.01GO  TO  2? 

IJJ-IJJ-21  • 

DO  23  I-l.IPN 
IJJSMJJ+I 

un.  i)  -HZ(i.u$)*pci»t) 

25  U(I.2)«HZ(IJJS) tP(I.2) 

GO  TO  29, 

2?  DO  28  1*1. I PN 

1JJS-IJJ+I 

U(I.n»2<IJJS)*P(I#l) 

20  U(  1*2)  *Z<  I.1JS)+P(I.2) 

29  DO  30  1*1.2 
1,11  *1.111+1  + 1 SET 

30  LRITE(3*  UlXLKJI.n.Jl-l.lPN) 

IJKM* IJ+l 

CALL  STASIK  IJKM. I rSPR) 

DO  34  K-1.2 
|-:i:-lk.iD<(lJK3-I0.2)+l 
KK 1 *3-K 

J J2*2+(KK- 1 > +KKI+ISET 
READ(3‘  JJ2)  aJCJJl.  D.JJl-1.  IPN) 

CALL  loop<iph.i:,y.w> 

74  CONTINUE 

UP  I TE  ^ 2 ) X 

1F(  1,1. LT. 3) GO  TO  36 
DO  35  I-l. IPN 

35  V<IJ»>«I>**2+V<!> 

36  DO  3"  I-l.IPN 

3?  y<  i l-o.O 

lFflC.LT.Fn.ONjf.i  TO  24 
40  lFMill.-nm  I.GC.IF1UG0  TO  50 

IF  I J.':n  i iPEADf.  l.tHD-r'51  (ZdP! . IP-22.  IPN221) 
On  aiUTIIIi.ii. 

r,n  to  oo 
re.  ip;,  rm  i 
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80 


100 

200 

300 

400 

500 

600 


END  FILE  2 
END  FILE  I 
END  FILE  3 
L|PITE(.5.500)V 


IFC  IPiitftT.EO.  I M,'F'1TE<5,600> 

S!;i;ms,e'  "**  °utp“f 

FOPirm*  ENTER  IFILV) 

FORMitK  10 15) 
fopiiiiTi.  i:;,  iciei2.4) 

FORMOK*  1END  OF  FILEV) 

EHt> 


DATA  FILEV) 


Figure  A. 18  (Continued) 
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DIMENSION  HC 1021 > 

CALL  LAG JNT(2«255»H) 

CALL  SGTF It U . ' LAGR . HAT' , IERR, ' DK' . 0> 
UR  1 TE  ( 1)  < H < I > , 1 » 1 , 5 1 1 ) 

END  FILE  1 

CALL  LAG  I NT <4, 235, H) 

CALL  SETF IL ( 1,' LAGD.DAT’ -IERR,' DK',0> 
WRITE* 1>H 
END  FILE  1 
END 


Figure  A . 19  SETLAG 
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C0m0N/DAT/X(255),0DAT<91 . 1? ) . PDATO  I , 

r mens  ion  isna.S).  ifi<2.5>,name(5) 

DIMENSION  FAKE (2040)  • 

DATA  IVAR.  ICSDR,HAMEU,',lS,75r0.3*'  U'.D'.'ATV 
URITE<6. 100) 

PEAP(G . 200)  (lint E ( 1 > . 1 - 1 ,3) 

WRITE <6. 30O) 

HEAD (6.400) 1BEG. IEND 

END  FILE  6 

IJK-0 

IF ( ( IEND-  I0EG)  .LE.  127)G0  TO  5 
JDIF-233+IBEG-lEND 
IIFIL-IDIF  + l 
IF(HFIL.GT.85)NFIL-85 

I TEMPS -IEND 
GO  TO  S 
IDIF-IENIt-IBEG 
NF1L-IDIF+1 

IF  tNF  IL .GT.85)NF IL-85 
ITEMPS- IPEG 

itempf-  i rcirs+16 

IJK-UK+1 

IF(  I J,K.GT.5)G0  TO  18 
IFUDIF.LT.  17)  1TEMPF-ITEMPS+ID1F 
IS1  ( 1 . 1 JK)  • ITEM’S 
I D IF  - ID  IF- 17 

IF ( ITEMPF .GT.255)  GO  TO  8 

IF  1 <1, 1JK>- ITEMPF 

IFUDIF.LT. 0)G0  TO  10 

ITEMPS- I TEMPF+1 

GO  TO  b 

IF  1 ( 1 < ] JK) -255 

ISK2.IJIO-1 

1F1(2» I JK) -1TEMPF-255 

IFUDIF.LT. 0)G0  TO  10 

ITEMPS- IF  1 (2.  IJIO  + l 

GO  TO  6 

IF ( I JK.GT .5) IJK -5 
WRITE (5,dG0) IS1. IF  1. IJK 
END  FILE  5 

CALL  SETFILU.* DELDOP.DAT'.  IERR,’DK'.0) 

DEFINE  FILE  l CNF IL. 120. U. I VAR) 

CALL  SETFILU, NAME,  lERR.'MT’  .0) 

DO  12  I - 1 -IIP  IL 

MRITEU  * 1VAR)  (X(J) , J-l  .64) 

DO  70  KL-1.1JK 

IND-M.L-n>l<t7 

1JITH.L 

CALL  STASIK I JLM,  ICSDR) 

C ALL  LOTA  i ' 1 . 1 7,  IS  1 ( 1 . KL ) . IF  1 ( 1 , KL ) . I ERR) 

I I **  1 

I."'-Uf  1U.I  LI-ISU  I.KD'M  l 

IFUSl'2.  :i.)  .HE. CO  12"  1 2-*  C IF  1 (2.KL)  - IS  1 (2, CL) ) + 1 
IFUEKP.IA.niCiA  TO  i" J 


Figure  A. 20  JDSPEC 
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lFfl.L.EO.  t)GO  TO  28 
GO  TO  22 

20  CALL  LOTATE*  19.64.  IS1  < l.KL) » IF  1 ( 1..KL) . IERR) 
IF ( 1ERR.EO.O)GO  TO  24 

22  DO  23  I-1EPP.64 
' DO  23  I J"  II . 12 

OOATH.I  D-0.0 

23  PI'AT(MJ)«O.0 
GO  TO  20 

24  DO  25  1-65,01 
DO  25  1J-I1.I2 
13“  I— 17 

CIDAT*  1. 1 J)  -QDATU3. 1 J) 

25  PDATt  I,  I J)  -PDAT<  13..I  J) 

28  SERP-FLOAT*  IERfi-D/64. 

IF<SERR.LT.0.)SERR»1. 

IFCSERR.EO.O.JGO  TO  70 
DO  50  I- 1 1.12 
DO  30  J-1,32 

QDATU.n-PTiATCJ.imJEJGHTtJ) 

30  PDATM,  I)»PDAT(J.  I)*UEIGHTCJ) 

DO  35  J-33.64 
J1-65-J 

QDATCJ,  D-CIDATHJ.  I )*UE IGHT(J 1 ) 

35  PDATCJ.  I)«FDAT(J.  I )*UEIGHT< J 1) 

CALL  FFT*QDAT( 1.1), PDAT* 1 , 1 ) . 64, 64, 64. 1 > 

DO  40  J-1,64 
TEMPI “QDATCJ. I) 

TEMP2-PDATCJ.I) 

40  QDATCJ,  I)  - ( TENF 1 #TEMP 1 +TEMP2#TEMP2) #SERR 

IVAR-IMD+I 

READ* 1* I VAR) CXCMP) .MP-1,64) 

DO  45  J - 1 *64 

45  >:<J)-X(J)+ODAT(J.I> 

I VAR ■ I VAR- 1 

50  I.,IPITI-<1‘ IVAR) CX(MP),MP«! -64) 

DO  60  IJ-1.  17 
DO  60  1-11.12 
I.U-lJ+64 

ODATt  IJ.  n-ODATClJl.  1) 

60  PI-HHIJ,  1 1 -PI'AT C I J 1 , I ) 

IFC  lEPP.GT.OGO  TO  70 
GO  TO  20 
70  REWIND  2 

END  FILE  2 

I VAR- 1 

DO  80  1-1,111- IL 

REAP'  1 ‘ I VAR''  CXCMP)  .IT- 1.64) 

i, IP  I TF  i 5.  r.nn  i i 

-jo  nr:  1 T|:; i r_i . vino . tip)  ,Ml'»l  .64) 

r.HD  ni.r.  5 

' EIID  TILE  1 

I CM)  » I.I|-,I  M r i ■ I'llTEF*  II:  lit*  OF  INPUT  DATA  FILEV) 
200  F0PllnT(3A2) 


Figure  A. 20  (Continued) 
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300  FOmvrr.  enter  imdicies  of  left  AND  right  top  boundaries*/ 

r MAXIMUM  CIRCULAR  SEPARATION  OF  TAPS-84'/) 

400  F0PIHK214) 

500  F0FIIATU0EI2.4) 

600  FORMATt-v'  RECORD  NUMBER-' 13) 

END 

SUPPOUTINE  LOTATEdl.  12-  IS t > IFJ.  IEPR1 

CPI 1!  mi  I 'DAT-'X ( 255 ) , ODAT < 8 1 . 1 7 ) . PPAT ( 8 1 . 1 ? ) > WE  IGHT< 32 ) 

DIMENSION  I3U2).IFI(2).I0FF(2) 

IERR-0 
IOFF  d ) "0 
N-2  . 

IOFF  <2)-(IFl<n-ISm>)+l 
IF< 151 (2) .EO.0)N»1 
DO  20  1-11.12 
F:EAD(2.END-5)X 
GO  TO  8 
5 I ERR  » 1 

GO  TO  25 

8 DO  lO  1J-1.N 

J5TART-IS1 ( I J) 

IFIH-IF1C1J) 

DO  10  J-ISTART. IF  IN 
J 1 ■ J- ISTART+1 + IOFF  < 1 J) 

10  QDAT d . J 1 ) *X< J) 

C END  OF  FILE  CANNOT  OCCUR  HERE  IF  IT  HAS  NOT  ALREADY  HAPPENED. 
READ (2)X 
DO  15  IJ-l.N 
1START-1S l ( I J) 

IFIN-IFldJ) 

DO  15  J-1START. IFIN 
J 1 -J-ISTART+l+IOFF ( I J) 

15  PDATd.Jl)— X(J> 

20  CONTINUE 

25  RETURN 

END 


BLOCi:  DATA 

COI II  ION  DAT (255) . ODAT (8 1 . ID  ,PDAT(8 1 . 17) .HEIGHT ( 52) 

DATA  X.PDilT,  PDAT.UE IGHT/3009*0.00,  . 2543.595.  . 2837055 . . 3 1390 16,  , .3448 

1242. . 37653631.. 4084023. . 4408209. . 4734933 . . 506290 1 . .53'?0r07, .571723 
29  ■ .6040391.  .6360360.  .6674292.  .6981291,  .770001 1,  .75621  1'7,  .7047307,  . 

33113316. . 8365931. .8603970. .0026353. .9032019. .9220010. .9389438. .953 

49502. . 9669485.. 9778773. .9966037. .9933259.. 99777 16. 1.0/ 

END 
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DIMENSION  SDATC64. 14)„MAME<5> 

CALL  SETFILd, 'DELDOP.DAT',  lERR.'DK'.B) 

WRITE <6, 100) 

READ(6,20fl>NFlL.NREC 

DEFINE  FILE  1 <NF IL> 128.U, I VAR) 

I VAR ■ 1 

NRECL-NFIL-NKEC+1 

NRECL-MAX0CNRECL.14) 

HREC-NREC-1 

IF (IIREC . EQ.OGO  TO  8 

DO  5 l«l,NREC 

5 READ ( 1 1 I VAR ) CSDATCJ, 1) . J«l,64> 

8 DO  10  l-MIRECL 

10  READC1’ I VAR) (SDATCJ, I)»J»1.64> 

DO  12  I -33,64 

12  LIRITEC5.30O)  CSDATC I , J) , J- 1 . NRECL) 

DO  13  1-1,32 

13  tJRlTG<5,30O)  (SDAT<  I . J) , J«  t » NRECL) 

100  FOPMATC  ENTER  NO.  OF  RECORDS  IN  "DGLDOP.DAT",'/ 
T AND  FIRST  RECORD  TO  BE  EXAMINED.'/) 

200  FORMAT <2  IS) 

300  F0RT1AT(  14E9.2) 

END 


Figure  A. 21  DDREaD 
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PIHEN3I0N  XOS) . D0PMARC64)  , PELIIAfi  (83)  .PELPK(S) . DOPFK  CIO). 
lIPELrK<5>.SPOPPKMB>.f>MiLD<5> 

DATA  POPI  IhR - DELI  WRv  1 49+0 . 0/ 

i jp  I TE  ( 6 . inrn 

PEAP(O.20n>HFIL 

Clip  FILE  9 

END  FILE  6 

CALL  SETF IL < I > ' DFLDOP.DAT' . IERR. ' PK' , 0> 

PFC I HE  FILE  KMFIL.  129. U,  1VAR)- 
IVAR-l 

HZPES’SOOO./'i  2S5.+64. ) 

HZQFFl— 32.+HZRES 
HZ0FF2  »3 1 . +HZRES 
PELRCO-l.^OOG. 

PELTOT-FLOAT  ( MF I L- 1 ) /9O00 . 

DO  30  I-l.HFIL 

REAPC 1 ' I VAR)  (X< J) . J-33.96) 

DO  5 J ■ 1 . 32 
5 X(J)«X<J+64> 

DO  10  .1-1.64 
SMfJL-l. 

IF(<J.EO.  lKOR.  0J.EQ.64))SMUL-.3 
10  DELMhR ( I ) ■DELI’IAR ( I ) +X(J)#SMUL 
SIHIL-1 . 

IF( U .EQ.0) .OR. < I .EO.MFIL) )SMUL-. 5 
DO  20  J - 1 .64 

20  DOPmRU)-DOmiR<J)+X(J>*SMUL 
30  CONTINUE 
END  FILE  .1 

UR  I TE ( 5 . 300 ) NF 1 1.  * DELRES 

LSI  I TE  ( 5 • 600 ) (PELI1ARU) » 1-I.NFIL) 

UPITE<5,500) 

UR  I TE  < 5 . 400 ) HZOFF 1 , HZ0FF2 . HZRES 

URITEUS.GOOiDOPMAR 

UR  I TE ( 5 . 50O ) 

DO  40  I-I.HF1L 
DEL  HIT-DEL  IimDELMARCn 
40  DELI  l-t'ELI  IlFLOATC  I)  iPELllARt  I) 

DELU-DELM/DELIHT 
DELI  15-PELNHI'ELRES 
DO  45  1-1.64 
D0PIHT-D0F'INT4D0PIWR<n 
45  DOPM-DOF‘1 1 FFLQATC  1-33)  +DOPNAR  < I ) 

DOPM-POPtl  POPINT 
D0PI1H-D0P1HHZRES 
DO  50  I-l.HFIL 

I F t DLL  MOM . I.  T . PFLHOP  < I ) ) DEI.MOX-PELUOR  C I ) 

LU  M'I.VHI;»1t:i.  V,  KlT.OilT  ( 1 1 -DELI  I J I- 1 2'1-DELI  IAP  (I ) 

DEL.  VHP  --DEL  '-VlR  I'FL  I M T 
DO  55  I -1.04 

t r u ion  io: : , i . i . Dornop  c n ) Don  in;  mdopmop  f i ) 

5,5  DiipyiiP»Ii0Pviirn-fr;L0rtTn-33)-D0PI1)-i-+2'l:D0PMAR(n 
pupviip  -Dorviir'  dop hit 

DI-I.5PS “SOP  ! U'ELVAR)  +DELRE5 


Figure  A. 22  DDSTAT 
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58 

eo 


65 
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65 
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75 


108 
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• trio 
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r,i.  in 


Dorr, dh-gmpt  < rmpvnK- > \ hjirfs 
up  m:  * s . r un  i deli  i'.. . m:i  sds 
MP 1 1IM  5,508 ) 

UP  I Tf£  f 5 . 800 1 DQPI  IH , DOP'aDH 
Mil  PIM  I ) "DEUINPC  |) 

I'OPPI  C 1 ) -1>0PI  INF!  ( 1 ) 

J«1 

I DEL  *0 

DO  60  1-2.IIFIL 
D I FF  -DELMNR  ( I ) -L'FLMAR  < I- 1) 

IFCDIfP.GT.O.)  go  TO  50 
1FCDIFF0.LT.0.)  GO  TO  58 
IF(  (DEU'fNRC  1-1  ^DELMAX)  .LT.  1 ,E-2)G0  TO  58 
DELPIMJ)  -DELMARC  1-1) 

IDEL-JDEL+1 
IDELRKC J) -1-1 
J • J+ 1 

DIFF3-DIFF 

CONTINUE 

J-l 

DIFFD-O. 

IDORr.0 

IDEU-IPEl.-l 
I F < l DELI ,60,0) GO  TO  66 
DO  65  1-1.IDEL1 

SDELD ( I > -FLOATC  1 DELPK ( 1+1) - IPELPKC  I) ) *DELRES 
DO  70  1-2.64 

D I FF -POPHrtR < 1) -D0PMAR ( I - 1 ) • 

IFCDIFF.GT.OM  GO  TO  68 
IF<D1FF0.LT. 0. ) GO  TO  68 

IFC  CDOPNARC  I- 1 )/POPIIAX)  .LT.  1 ,E~2)GQ  TO  68  . • 

poppi:i:j>»DOPmR<i-i> 

IDClP-IDOP+l 

SDOPPIM  J)  “Fi.OflTC  1-34)  *HHRES 
J-J  + l 

D1FF0-D1FF 
CONTINUE 
WRITE (5. 500) 

UR  I TEC  5. 90C0 

WRITE rg. 1000) ( IDELPKC I) . 1- 1 . IDEL) 

IF  C 1I>I:L1.EC!.0)GO  TO  75 
WPITC5.  1200) 

WRITE (5. 600) CSPELPCI), I-l. I DELI) 

UP  I TEC  5. 580.' 

WRITE  C5. 1 100) 

WRITE C5.S00)  CSPOPPKCI),  1-1.  IDOP) 

FORlKlT ( ' ENTER  HO.  OF  DELAY  TAPS  PROCESSED*/) 

FOPMilTf  T*M 

rori r. iTi- • i ii.mc.iiihi  dfi..hy  power  disipidlition  fur  tiif  -.it, 
iMIHFD.  C iiC  M TOP  18  SEPARATED  PV.EI0.3.-  KECnim*!.  V) 
roi-r 1. 1 rr  • umi'i.iiimi  duppuip  puufp  pisipiuut imn  i-ptiim  .r-rs  • 
iff.?.’  1 1,*'  in  sups  fjr * . r r . 3 . * irvi 

•F'OITIn  I''/.  . i 

rrjpiinr.'j::.  i nr.  i . t:i 
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* Tors  n::n 

• n;c  in  j 


70  0 


FORWTt*  DELAY  1CAH«'E10.3*  SECONDS COPRESPOMD IMG  STANDARD  DEVIA 

1TIOII='EI0.3*  SECONDS'/) 

SOO  FORI  In  * f ' DOPPLER  1EAII»*FI0.3'  NZ CORPESfOtlPlNG  STANDARD  PFVIA 

1TI0H»'F10.3*  II Z’S) 

SIH)  FOPI  V|T ' DELAY  PIAKS  OCCUR  AT  TOPS:'/) 

1000  FOPDYn'jX,  101S> 

1100  FURirm*  DOPPLER  PEAKS  OCCUR  AT  FREOUEIICIESUN  112):*/) 

1200  FOPmTC  HIT! I CORPESPOtlDING  Tilt  DIFFERENCES  111  SECOND?  RETUGEN  I»C 
1LAY  PEAI'S  OF:*/) 

END 
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DIMENSION  Y(0??O) . VS1 ( 12S) . V3P(255) . VS2< 123) > VSPC255) .POKE (2048) 
EOU I VALENCE !<  YU ) - VS  1 ( 1 ) ) . < VS  1(1 ) . VSCI  ( I) ) . ( VS2  ( 1 ) . VSP  U ) ) 
C0MM0N/PANSID.-I1.  12 
1 1 - 

1 2 —243  "6 
WRITE  <6. 100) 

END  FILE  6 
PEAI’ (0.300)  I F 11- 
END  FILE  8 

CALL  SETF IL<3.  'HFS1MS.DAT' . lERR.'MT' .0) 

CALL  SCHAM(VS1.VS2) 

1CSDR-" 177370 
1FIL-1FIL+4 
DO  50  J-1.IF1L 
IJKM-J 

CALL  STASIK  IJKM. ICSDR) 

CALL  CHANCVS1.VS2) 

VS1CD-0. 

VS2CD-0. 

CALL  FFT(VS0(1),VSP<1). 233. 253. 253.-1) 

DO  10  1-1,255  ■ 

10  VSCK  I) -VSCK  D/255. 

CALL  FFT(Y( 1 > . Y<2> . 384. 384. 384. 2) 

CALL  REALTR(Ya).Y<2>. 384.2) 

URITE(3)Y 
DO  20  1-256.770 

20  Y<n-0.0 

DO  30  1-129.255 
Y(  I) -0.0 

30  VSPCD-0.0 

50  CONTINUE 

END  FILE  3 

100  FORHiVTY*  ENTER  IFILV) 

2O0  FORHATUX,  10E12.4) 

3O0  FORMATt 15) 

END 

SUBROUTINE  SCHAIKVS1.VS2) 

COMIION-'RAUSID/1 1, 12 

DIMENSION  Y(?f .P(4),B(4).FPQP(4).TD<4) 

DIMENSION  VS  1 ( 128) , VS2( 128) 

REAL  K0 

COItl  IQIN'F  ILT-'A  I (.4)  • A2  (4) . A3(4) 

COftlON/CHAI  ID-NPATH . NTONE 

COMMON, CHAIN',  CHPH(4) .FIXPHC 128.4) .X0(8).X1 <B) 

DATA  Tl JHP I /S . 283 1 853'' 

WRITE  05.  109) 

RE AD f 8, 1 12)  IITONE 
WRITE -.6. 110) 

PEAD'A.  103)FSEP 
IP  I IE  ft. 105) 

PCulHO.  103)  PATE 
. LIP  I IE  ft..  100) 

READ  . 1 12)IIPhTH 
DO  5 l-l.NPATH 
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writer. 102^  r 

READ'S, 103)B(J) 

WRITE'S, 104) I 
READ'S. 103)P(1) 

WRITE (6. 107) I 
READ'S. 103) FM)PC I) 

WRITE'S,  IIP  I 
PENDU3.  103)  men 
5 CONTINUE 
OFFSET"0. 

TFRAI  € « 1 . /ROTE 

DO  15  M.MPATH 

PU)«10.**C(PU)-3.01)/I0.) 

KO -EXP  ( -6 . 233 1 853*8 (I ) *TFR AME ) 

A1  ( I) »2.*K0 
KO*K0#*2 
A2'  I)— K0 

A3' I)  -SOPT'P ( I ) n l . -k'0)**3/<  I . +K0 ) > 

RHO-Al(I)/(l.+K0> 

1*0  15  J-1.2 

1F(ABS(B( I ) )- 1 ,0E“6)B.8,9 

8 V< 1 > *1 . 

y<2)«l. 

GO  TO  10' 

9 CALI.  GAUSZ(Y) 

10  X0(2*!+,t-2)  "SORT<P<  I)  )*Y(  1) 

XI  (2  KI+J-2)  ■SORT(P(  I)  )*(RHO*Y(  1 >+SORT(  I .-RH0**2)si<Y(2> ) 

15  CONTINUE 

DO  16  J-l.NPATH 

CHPH<  J)  -TLIOPI»FDOP(J)*TFRAME 

A«TU0P1*FD0P(J)*TD(J) 

DO  16  M.NTONE 

C-TWOPI*TD'  J)  #(  (FLOAT'  I)-l.)  *FSEP+OFFSET) 

F I XPH ' I , J ) -HMOD ( A+C , TUOP I) 

16  CONTINUE 

CALL  CHAN < VS  1 »VS2> 

RETURN 

100  FORMAT' ' ENTER  NUNBER  OF  PATHS  (FOUR  OR  FEWER)'/) 

102  FORMAT' ' ENTER  RMS  BANDWIDTH  OF  PA3HM2.'  IN  HZ'l)  . 

103  FPRMAT'F 14.9) 

104  FORMAT''  ENTER  POWER  OF  PATH’. 12. * IN  DB'/) 

105  FORMAT'’  ENTER  FRAME  RATE  IN  FRANES.'SEC*/) 

107  FORMAT1;*  ENTER  DOPPLER  OFFSET  OF  PATH*.  12.’  IN  HZ'/) 

108  FGPMAT ( ' RAN  TERMINAL  COUNT' . 2 ( 18) ) 

109  FORMAT' ' ENTER  NUMBER  OF  ACTIVE  TONES'/) 

110  FORMAT''  ENTER  TONE  SEPARATION  IN  HZ'/) 

111  FORMAT'’  ENTER  DELAY  OF  PATH'. 12.'  IN  SEC'/) 

112  FORMAT • 15) 

END 

'■•  UDROUTINE  nil'll  1 1;  VS  1 . VS 2 ) 

Mill  1 1- . I PI  I V"M<  i::B»,VS2'123) 

, c oi  ii  ii  h i.  phi  is  inn,  12 

MIK'H'iHill  Vi .?'•  .P(4«  .AiTIPiOP 

an  men  1 nimiD  iipath.ntone 
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cor  It  ion.  C HmIIIi  -c  MPH  C4>  ,F1>;PH<  128,4)  .X0(8),X1(B> 
DATA  Tl.WPI/15. 2831853/ 

CALL  F1LT2(IIPATH,X0,X1) 

DO  10  J • 1 » NPATH 
K-2YJ-1 

R ( J ) -SORT < XO ( K ) **2+X0 < K+ 1 ) *#2 ) 

G Al  IMA  < J ) • ATAl  12  C XO  ( K+ 1 ) , X0  ( K ) ) 

10  CONTINUE 

I'fl  20  1-1,  HTOHE 
VS  1 C 1 ) -0 . 

VS2(J)-G. 

DO  15  J-l. NPATH 

F I XPH  (I . J ) -PI  lOD  (F t XPH  < l , J ) +CHPH  (J),  TUOP I ) 

A-F IXPH( I . J)+GAI1MA(  J) 
VSKD-RCJMCOSW+VSKl) 

VS2( I) -R(J)*S !N(A)+VS2( I) 

15  CONTINUE 

20  CONTINUE 

RETURN 
fcMD 

SUBROUTINE  F ILT2(N.X0.Xn 

COmWlVRAHSID/Il.  12 

COMMON  ■•F1LT/A1  <4)  ,A2<4)  .A3<4) 

DIIENS10II  J‘SH8).Xi<0).Y<2) 

DO  10  1-1, N 
CALL  GAUSZ(Y) 

DO  10  J-l, 2 
K-2*l+J-2 

VC J) -A  1 ( I ) *X0 <K) +A2 C I ) *X1 (K) +A3 ( 1 > *Y<  J) 

XI OO  -X0<K> 

xoao-vu) 

10  CONTINUE 
RETURN 
END 

SUBROUTINE  GAUSS (Z) 

COMMON-PAHS  ID-'I 1. 12 
Dll  (FUSION  Z(2) 

x-piiic:  1 1 » 12) 

Y-rUIIU  I,  12) 

;>sor'T(-2.iAi.OGi:x>) 

Y-S.  I'OSIGSS-HY 
2(1' -XICOS(Y) 

z(2)-;msiiky) 

RETURN 

END 
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intrusion  Y'rrg)  ,2(765) .LH770) .H(51 1> ,P(5I0) .X1255) 
DIMENSION  FAKE (2048) 

EPIJ I VALENCE  <H(2).P<1>) 

COMPLEX  YCC303) ,UC<385) 

EQUIVALENCE  CYCn.YCil)),  (UCn.UCU)) 

DATA  Y/770i;0.OX,  2/765*0. 0X 

DATA  IQ.  II.IPII.  ICSDR/2, -1,255, "177570/ 

DATA  1PN2, IP1I3. IPN4/5 10.765. LO20/> 

SIPND4-1 ./(4.0*FLOATC  IPN+1>*768.0> 

ISET-0 

CALL  SETF IL (1 , ' LASR . DAT' , IERR. ' DK ' . 0) 

READIDH 
END  FILE  1 
WRITE <6. 300) 

END  FILE  6 
READ! 3. 100) IFIL 
END  FILE  0 

CALL  SETF IL< 1. 'PNSEQT.DAT' . IERR. 'DK' .0) 

DEFINE  FILE  1 (1. 1540. U. IVAR) 

CALL  SETFILO.  'HFS1MS.DAT' . IERR.  'MT'  ,0) 

CULL  SE  IK  I L «: 3, ' TESTD . DAT' . IERR . ' DK'  .0) 

JFIL-IFU.+4 

DO  70  IK-1. IFIL. 2 

DO  70  IXY-1.2 

I JKM- IK+IXY- I 

CALL  STASH!  IJK11.  ICSDR) 

READ ( l ’ 1)  Y 
PEAD<2.EN.0-75)U 
DO  30  1-1.395 
30  •it<I)-YC<I)*LC<n 

CALL  REALTRCY< 1) ,Y(2) .384. "2), 

CALL  F*T(Y( l ) . Y(2) .384.364.384,-2) 

? !-l!0D(<Il  + l).3> 

1 1 1 ■ 1 1 fcIPN 

IF  (Ill.QT.IPM)  GO  TO  34 
DO  32  I-1.IPH2 

32  Zi l + Iin-Yri  + IPN)*P(l)+Z(I  + ni) 

GO  TO  38 

34  1 14« IPH3- 1 1 1 
DO  35  I-1.II4 

35  El  ; I + t 1 1 J-Yd  + IPN)*P<I)  +2(1  + 1 ID 
1 12  - IPN4- 1 1 1 

113-  IPII3-  1 1 1 
1.14*1 1 1-IPM 
DO  37  !»!.  1 14 

37  Z>.  I • - I + II?.m.P(I  + II3)+E(1> 

38  DO  39  I-l.IPN 
IMH111 

39  .7*  I l<l)  *,"’i  I I P ISIPMD4 
If  ■ If  .1.  I . ? u,i.i  If.i  .10 

ir  i i::r.iiL.  t • go  to  43 

DO  .11  I»i . H-'tl 
4i  > : • i 1 * 1 1 + 1 1 1 ) 

GO  TO  45 
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43  URITEC3) (X'J).J-l. 1PN).  <Z< 1+111 >. I»l, IPN) 

45  DO  50  1-1. IPN 

50  Zd+im-0.0 

70  CONTINUE 

GO  TO  80 
75  1 SET" 1 

80  END  FILE  3' 

El II'  FILE  2 

END  FILE  1 

1FUSET.EG.  mJRITE(5,500)IJKM 
100  .FORM0TU5) 

300  FOPMhTt'  ENTER  1FIL'/) 

500  FORI IrtT ( ' END  OF  FILE  ON  RECORD' . 15/) 

END 
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DIMENSION  YI770) 

DATA  ISTATE,  ITAP/'l*  "270^ 

Y(l)«FLOAT(  ISTATE) 

DO  15  1-2,255 
ISTATE*  I FSF!<  ISTATE.  1TAP) 

1FC ISTATE. AND. 1)5,5*  10 
ITEIt5-- 1 

Y( I ) -FLOATC l TEMP) 

I TEMP -1 
DO  20  1-1,255 
Y<  1+255) -Yd) 

YU+510)-Y(t) 

URITE(5.500)Y 

CALL  FFT(Y.Y(2). 384.384. 384,2) 

CALL  REALTR(Y,Y<2) .384.2) 

CALI..  SETFILU.'PNSEQT.PAT*.  IERR.*»K'.8) 
DEFIII'.:  FILE  1 Cl,  1540,1).  IVAR) 

UP  I lECl'DY 
END  FILE  1 

CALL  REALTRCY.Y<2). 384.-2) 

CALL  FFT(Y.Yt2). 384,384,384.-2) 

URITEC5,5O0)Y 

FORI-MTC IX,  10E12.4) 

END 


1 
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S'JErOUTJNE  FFTCA.B.IIIOT.N.NSPAN.  tSN) 

mifikion  <\< i ) ,D' n 

MIEIIOIOII  lll-TiCC'*.-* 

DIHF.H8H.1II  ftT(5)»BT(Sl 
FOMIvVlLEtlCE  (1. 1 1) 
tH',F-5 
MA>F«23 
1-0 
JF-0 
K-l 
11*5 
KT-l 

MFflC ( n -4 
M(mC<:2)-4 
HFf»C<3)  ”3 
HFhC<4)-2 
HFAC',5)  «4 
INC- ISII 

PAP«8.0+ATAN<1.0) 

S72-FAP/5.0 

C?2-COS(S?2) 

Sr2-3IH(&72) 

S 12Q"SORT(0.?51 

IF  ftSII.GE.O)  GO  TO  10 

S72— S72 

S120-SI20 

PAD-RAD 

IHO-IMC 

10  MT*1HC*HT0T 
kS-IHiMHSPAN 
KSPAH-KS 
Htl-MT-IHC 

jc.m.s/m 

RADF-RHp*FLOnT(JC>#.5 
100  SD-F'liDF  'FLDAT(KSPAM) 

C D " 2 . 0 +S I M < G P ) ,t!*2 
SD'S  IlKSD+SD) 
i:r- 1 

i-i+i 

IF  OIF  AC.  I)  .HE. 2)  GO  TO  400 
NSPAII-KSPAII'S 
K.  1 »K8PAH+2 
210  K2*I':+KSPAN 

hi;-0(K2) 

BK«B(K2> 

A<K2>-A<KK)-AK 

EKl;2)-B(KK)-BK 

H <1.1  i -AOOO+AK 

BO. I > 'BO  t )+BK 

I K-I\2-»KSPAH 

IF O I .LE.Itll)  GO  TO  210 

M 

IF  OI.LE.JC)  GO  TO  210 
IF  O I .GT.I  SF'GHj  GO  TO  GOO 


Figure  A. 26  Subroutine  FFT384 
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r 

r 

jC 

L 
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220  CI-I.0-C1' 

SI -SI' 

230  I 2*1  T »kSPAN 

Ah  •-rtO'2'1 
bi:-D(i:K)-D<K2) 
Afi:n«iUm+A(K2) 
B'KI:>-Bi.KIO+B(K2) 
A(h.2>*Cl*iU'-Sl  '«K 
Et'  i:2  ' -SI+hIC+CH-BK 
IT  *k  2+h  SPAM 

ir  tl.T.LT.HTl  GO  TCI  230 

l'2-KK-HT 

Cl— Cl 

KK-K1-IC2 

IF  (klC.GT.K2)  GO  TO  230 
AK  »C 1 - ( CMC  1 +SL*S 1 ) 

S 1 * ( SMC  1 -CM'S  i > +S  1 . 

Cl-AK 

ir 

IF  (KK.LT.K2)  GO  TO  230 
K1-K1+1I1C+INC 
KK»  (1C  1-KSPAN)  /2+JC 
IF  (KIC.LE. JC+JCI)  GO  TO  220 
GO  TO  100 
320  K 1 ■KK+KSPAN 

IC2-ia+l(SPAH 
AK-ACKK) 

BIC-B'IT) 

AJ*h(IC1)+A(K2) 

BJ*B<K1)+B(K2) 

A(KIO  »AK+A,I 

BOCIO-BK+BJ 

Al'  — .5-tAJ+AK 

DK— .5l«J+BIC 

A>CAnCl)-A<IC2))*S120 

6WB(M)-BiK2))*S120 

A(i:i>*Al>BJ 

A' 1.2)  'AK+BJ 
E!''K2)  ■£>!>- A J 
I I "1-2+ICSPAII 

IF  i I.I  .LT.NM)  GO  10  320 
k»k'K-HH 

IF  (I'K.LE.KSPAN)  GO  TO  320 
GO  TO  TOO 

400  IF  fllFAC'  I)  .HE. 4)  GO  TO  600 
I.SPIIII-KSPAII 
ISPAIl-kSPAU/4 
410  Cl'1.0 

01*0 

.120  1C  1*1  T-HTTAII 

12*1  1+1  SPAN 
h 4 J I I -c.l  •ill  I 
AI;P-A'Ki:)+H<K2) 
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hi  i r-n  < ii  * • m ■ i ‘.'.i 
h ir  - 1 1 ■ i i > ih'I  ?.  i 
H.m-ii-1  r>  m 
III  I l.'l  -Al.P+flJP 
AJPmiI.P-AJP 
Bl  P 'B'TJ.)+B<K2) 

eh 

tnF-nn:n  ibii<3> 

ruii-pii  i >•  -B'l'S) 

Ell M-.  I'Lil.P-iBJP 
Ei.lP'Eil.P-P  JP 

ir  n$n.LT.o)  go  to  456 

Mi:p»rtKH-DJn 

Pit  i i-nu  t*-ti  jm 

BKP-BKIH-AJII 

BI.TI-BKM-AJM 

IF  CS1.E0.O.O)  GO  TO  460 
iKKl.i-PKnCl-BFP^Sl 
B U<  1 ' AKP  +S 1 +BKP*C  1 
AU'.2)«AJP.iC2-BJP*S2 
ni'l-:2)-BJP:*C2H'BJPit€2 

ti  i K 3 • »ril  :i  HS3+BkH*C3 
M' ■U'3+KSPAH 

tF  EKK.LE.IIT)  GO  TO  420 
C2-Cl-(CDHC1+SD*S1) 

S 1 • ( SD*C  1 -CD'KS  1 ) +S 1 
CI-C2 

C2  "C  |ss+2-S 

$2»2.0>i>Cl>*Sl 

C3*C2i€l-S2ti51 

S3«C2l'Sl+S2*Cl 

IOKK-HT+JC 

IF  (UK.LE.KSPAN)  GO  TO  420 
KOKK-KSPflM+IHC 
IF  (KK.LE.JO  GO  TO  410 
IF  (I'.SPAM.EQ.JC)  GO  TO  800 
GO  TO  100 
AI.P -AKIH-B  Jl! 

AI.II-ALN-BJH 

EIP-tim-flJM 

BKII-Bm+AJM 

IF  CM  .HE .0.0)  GO  TO  430 

lV.K.n»AKP 

B(K1)»BKP 

A<T2  i «A  TP 

BH  2'»BJP 

A<l  3 ' -Hi  lt 

BU;3)-BI-:i1 

I i -I  i l ?FhII 

]!•  M.L.LE.IIT)  GO  TO  420 
GO  TO  440 
1 iirYici; 

I GPIIII-KSPAH 
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I M ull  I '.Pull,  i: 
ii  'i  ,i  ii, :)  r,n  in  T.’O 
if  1 1 .r.o.tti  CjO  io  boo 
1 1'" .11+ 1 
C2-I  .O-CP 

51- S1J 
C 1 -C2 
;-2>si 

|;|>|.|. +ISPAH 
Ml  *H  (KIO 

H*  1 1 > "C2  n'il'“S2>*D<KK) 
b 1 1 r 1 ioi  +c^'ifhi:im 
I ,l>  I I •ll  .'iPHII 

IF  • PI'.LE.im  GO  TO  730 
0K-SMS2 

52- $MC2+C1*S2 
C2-C1HC2-AK 
kk-IT.-IIT+KSPBII 

IF  (Kk.LF.KSPHM)  GO  TO  730 
C 2 -C 1 - ( ClM'C  i +SD*S  1 ) 

S1-S1  + (SD+C1-CD*SD 

i:k-i:i.-K3PMN+jc 

IF  CKIM.E.I'SPflN)  GO  TO  720 

I.OII-lsSPfiH+JC+IHC 

IF  (Kk.LE. JC+JC)  GO  TO  710 

GO  TO  100 

HPCD’KS 

IF  O'. T .EQ.0)  GO  TO  830 

i.-kt+i:t+i 

IF  (M.LT.IO  K-K-l 
J-l 

IIPlK+D-JC 

HPlJ  + l)»IIPlJ)/HFPC(J) 

IIPOO  -IIP(K+1)4NFAC<J) 

J-J+l 
K-k- 1 

IF  CJ.LT.K)  GO  TO  810 

P3-MP(K+1) 

l'.SF>‘ill"IIP<2) 

IJM.1C  + 1 

1',2-IGPrtH+l 

J-l 

Ml  «m  Cl'IM 
HOI-  .'-B(K2) 

H K2'  -BK 
OIMGO'.IO 
po:n»D(K2) 

BO'2)«8k 
I.OKIs  + IIIC 
l'2-l.?'PflU+l<2 

IF  0-2.LT. I- S'  GO  TO  820 
I 30  2-IIPU? 

J-J  + l 

1.2'- 1 ip  <:j  1 1>  +i-2 
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IF  ' l ’.nT.HP(J)  :•  GO  TO  030 
.1-1 

041)  IF  (I  I .1  T.1.2)  GO  TO  020 

i i>mmt  ■ 

K2»i  sr i'ih+i:2 

IF  (I  2.LT. I S)  GO  TO  e4f) 

IF  O.C.LT.KS)  GO  TO  830 
JC-13 

890  IF(2 iKT+1 .GE.M)  RETURN 
I. SPIII H IP  (KT+l) 

.1  "11-KT 
HFACt J+i'-l 

900  HFftC  U)  -HFftC  < J ) *HF  AC  < J+l ) 

J-J-l 

IF  < J.HE.KT)  GO  TO  900 
KT-HT+l 

uti*  nr  nc  m<t)  - i 

IF  fllN.GT.IIAXP)  GO  TO  990 

JJ-0 

J-0 

GO  TO  906 
902  JJ-JJ-K2 
K2-KK 
K-K+l 

kk •HFftC (K) 

904  JJ-KK+JJ,1 

IF  ( JJ.GE.K2)  GO  TO  902 
IIPC  J.’* -JJ 

905  K2 -HFftC (KT) 

K-KT+I 

KK -HFftC  (K) 

J-J+l 

IF  (J.LE.HN)  GO  TO  904 
J-0 

GO  TO  914 
910  P.-KK 

KK-NPOO 
HP(K)  — KK 

IF  (KK.IIE.J)  GO  TO  910 
K3-KK 
914  J-J+l 

M>HP<  J) 

IF  (IK.LT.O)  GO  TO  914 
IF  (Kk.HE.J)  GO  TO  910 
MP<  J)—  J 

IF  (J. HE. Nil)  GO  TO  914 
iih:  :f  "lift  xf+ me 

GO  TO  950 
924  J-J-l 

IF  (HP(J).LT.O)  GO  TO  924 
JJ.JC 
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II  . I J-l  SPAN 


Figure  A. 26  (Continued) 


923 


932 


936 


940 


950 


998 


MIPfJ) 

j I *jcU ,+11  + JJ 
i;  i »|  (.•tK'jPrtH 
12-0 

I 2*12+1 

ftTU2)»8<Kl> 

BT(K2) »B(K!) 

II  .|  l-l  HO  g>Q 

tr  (J-'.l.liE.KK)GO  TO  9.-8 

y t *|  K+K3P8H 
12-1,  t- JC*«.K+HP  f-K>  5 
l --IIP’OO 
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'irSIw-**'  00  ™ 036 

00  TO  932 

KM;K+KSPAH 
K2"0 
K2-K2+1 
ft(Kl)  «ftT(K2) 

EfKl) ■BT(K2) 

rF 

J-K3+1 
,|T»HT-K5PMH 

G0  T0  924 

return 

1SH-0 
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PEftL  HI 
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KSwT«M<T.O)^0BT(M) 
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5Ii-SUUSt,+SI,) 
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C ff'PL AV3  ericr  GEHEPATED  SNAPSHOTS  of  hf  channel  for  plvscf;this  is 

C *•  THE  OVERLAY  VERSION 

conmii  ovr- vorroi ..’(07bS) ,LKorro'.H(05i n .;:(25S) 

com  mu  • iMivni.  i i.  icsi»p.  io.  ipn.  iphiq.  iphioi  . iphic>2.  iset,  iopi.  ik,  iky 

M I CHS  I OH  P*0510> 

Mints  mil  FHI  E (2049) 

EQUIVALENCE 'H<2).P<1>> 

IPMIQ»1P» IPH 
IPHIOi  - IPIIIO+IPII 
IPHIQ2" IPHIOI+IPN 
I OP  1 ■ 10+1  . 

S 1PHD4-  I . 0/(4. 0i*FLOAT<  1PN+I  ) *768. ) 

IF IL- 10000 

CALL  LOAD ( 'LAG' .1, IEP) 

CALL  SETFIL  < 1 » 'PM5E0T.DAT' , IERR, 'DK' »0> 

DEFINE  FILE  1 C 1 . 1540. U. 1VAR) 

CALL  SETFIL  <2. 'HFSJIIS.DATMEPR.  MIT' .0) 

CALL  SETFIL  (3. 'TESTD.DAT'. IERR, *DK'.0> 

CALL  LWIT 

IF(IER.EO.O)GO  TO  20 
I JKH»* 100000 
GO  TO  80 

20  CALL  LAGINTC IQ. IPN.H) 

I'O  70  n;-l.IFIL.2 
DO  70  I XV- 1.2 
1JKN-IK+IXY-1 
CALL  STASH ( IJKI1,  ICSDR) 

IVAR-1 

CALL  LOAD (' SUB  1 0. IER) 

IF(IER.EO.0)GO  TO  25 

!SET-“2000O 

GO  TO  30 

25  CALL  SUB  1(1 VAR) 

IF ( ISET. HE. 0)GO  TO  00 
CALL  LOAM 'FFTOVR'.O,  IER) 

IF(  lER.EO.OGO  TO  30 

ISET»"4Q000 

GO  TO  80 

30  CALL  FFTOVR 

CALL  LOAD ('SUB2',0, IER) 

IF ( !ER.EO.0)GO  TO  35 
ISET- "60000 
GO  TO  80 

35  CALL  SUBS <3 IPND4) 

70  COHTIHUE 

60  El  ID  FILE  3 

EHP  FILE  2 
END  FILE  1 
I JKII-  IJI'I'H-I  SET 
CALL  STASH*:  I JKH.  ICSDR) 

EHD 

BLOCT  DATA 

COI II  WII  flVP-.-f  0770).. Tf  0763>,  l,K0770l.H(05in.)«255) 

COIIMUILDAVAL  II.  ICSDR.  10.  IPH.  1PII10.  IF'tIKlI . Il'llin.’.  ISET,  tPPl,  IK.  K.V 
, DATA  V.. -V  1535=1 0.0/ 

DATA  II. ICSDR, 10. IPH. ISET/- 1 . * 177570 • 2 . 255 • 0/ 

EHD 
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I.  I'll  lYni.r  in  i 1 ll<  -P(  I) ) 

1 1 im'm  Miin.  inpi  1 

ii  1 - in  i r 11 

IF  1.  Ill  .1-1.  I PI  I)  GO  10  34 
i»o  30  l-i.iPMio 

zu+nt  i»V'.n]rii)+p<i)+z(i+iii) 

GO  TO  30 
I 14*1 1 I'll  1 0 1 - III 
I'd  3*1  ni.ni 

.mu  1 1 »«vi  ii iPiinpi n+2<  i+i  in 
1 i:-ipiiio2-ii  1 

IU-IITII01  -III 
114  III  IPII 
t'f.i  MP  I -1.114 

t.1-,'1  HII2>  MJ(  t + 113>+Z(I) 
r-o  m 1 1 1 . ipii 

I I . 1-1+1 1 1 

2(1 14>  «i!i  1 14)  ISIPHI'4 
IF  1.  IK. LT.  I OP  I ) GO  TO  45 

1 r 1 i::y.!ie. nno  to  43 

I'G  41  1-1,  IPII 

•o  n-pt  nun 

GO  Til  45 

WP  HE  (3 1 CX( J) » J- 1 , IPII) , (2  ( I + 1 1 1 ) » I - 1 • I PM) 

I’O  50  I -1,1PM 

zdiiin-a.'n 

RETURN 

EIII) 
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' DIMENSION  X(5 12) * Y(512) .NAME(5) 

I DATA  NAME/3#'  '.'.D'/flTV 

URITE(6»300) 

IREnr.K6.400)  CMAMECI).  1-1.3) 

UR  ITU (6. 2O0) 

READC8.  100)  1REC 
END  FILE  0 

9 END  FILE  6 

I CALL  SETF IL( l.NAME. IERR.  'fIT'  .0) 

DO  10  I-l.IREC 
REAtK  1) (X(J) ,J- 1.255) 

I 10  READC1HY(J).J-1.255) 

| END  FILE  1 

MR  I TE (5. 500) (X(J) .J- 1.255) 

, URITE(5.500) CY<J). J-1.255) 

DO  15  1-1.255 

I 15  Y(D— Y(I) 

COLL  FFT(X. Y.512.512.512. 1) 
i DO  20  1-1.512 

20  X(I)KX(im2+Y(I)«)#*.S 

CALL  SETF IL < 1 . 'FFT255. DOT' . I ERR. ' DK' .0) 

UP  I TE  ( 1 ) X 

I END  FILE  1 

I URITEC5.500)X 

100  F0RMATU5) 

, 200  FORMAT C*  ENTER  RECORD  NO.'/) 

I 300  FORMAT ('  ENTER  NAME  OF  INPUT  DATA  FILE'/) 

, 400  F0RMAT(3A2) 

500  FCRMATC IX. 10E12.4) 

. END 


DIMENSION  X<513).Y<513) 

1 DATA  XMIM.XHrtX.YMIN.YMAX/'O..  1 . ,-60. » 10.>* 

DATA  XSIZE.YSIZE.'G..7./ 

CALI.  SETF IL  <1 . ' FFT255 . DAT' , IERR. ' DK* . 0) 

PEAD(1)(X<1), 1-1.512) 

EtID  FILE  1 
SHAX-0. 

DO  5 1-1.512 

X(  I)  -20. DIALOG  10 (X(  I) ) 

3 IF<X<I) .GT.SHAX)SMAX-X(I) 

DO  7 1-1,513 

Y( I) -FLOAT! l-li/512. 

X(  I)  -X(  I )-SI1AX 

7 tF<X!l).LT.~60.)X(l>— 60. 

X(513)-X<1) 

DX-  (XMAX-XI1IN)/XSIcE 
D Y-  ( Yl  IAX-YM I N ) /YS IZE 
XARG-XMIN/DX 
VAPlj-YlllN/DY 

CALL  MODE (2. XS IZE. 0. .1.5) 

CALL  MODE ( a.  YS1ZE.0. . 1.5) 

CALL  MODE <8.  >31  IN, DX. XARG) 

CALL  MQDE(9, YMIN.DY. YARG) 

CALL  DR AU ( Y ! I ) . X ( 1 ) . 5 1 3 « 44 1 ) 

XS12E1-XSIZE+.5 
YS1ZE1-YSIZE+.5 
CALL  MQDE!2,XSIZEl.-,75» 1.5) 

CALL  M0DE<3. YSIZE1.-.75. 1 .5) 

CALL  I'IODE(8,XM1N,DX,0. ) 

CALL  MODEO,  YMIN.DY, 0, ) 

CALL  MODE ( 7 . XS I ZE . YS I ZE , 9999 . ) 

CALL  AXES (37.1,' FREOUENCY  NORMALIZED  TO  SAMPLING  RATE' .23. 1 . 'MAGNI 
1TI.IDE  RESPONSE  <D9)') 

CALL  MODE < 8 . XII I N . DX . XS 1 ZE ) 

CALL  MODE  ! 9 , YI1 1 1 1 , DY,  YS  1 ZE ) 

CALL  AXES (0 . 1 , 'X* .0. 1,  * Y* ) 

CALL  DRAUKO..O.. 1.9000) 

cull  vers (0,0) 

END 
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investigator  has  some  knowledge  of  the  channel  conditions. 
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reconstruction;  (2)  prefiltering  to  improve  delay  resolution; 
(3)  channel  snapshot  generation;  and  (4)  Doppler  estimation. 

A set  of  flexible  Fortran  programs  which  meet  these  specifics* 
tions  are  described  in  detail.  Software  verification  is 
achieved  by  means  of  a program  generated  single  sideband  HF 
test  channel.  In  addition,  programming  changes  to  the  pre* 
viously  reported  channel  measurement  and  reproduction  software 
are  documented.  These  result  in  significant  decreases  in 
computation  time. 
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